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. Abstract 

O 

00 Highly spinning classical strings on M x 5 are described by the Landau- 

Lifshitz model or equivalently by the Heisenberg ferromagnet in the 

thermodynamic limit. The spectrum of this model can be given in 



/\ terms of spectral curves. However, it is a priori not clear whether any 

^ given admissible spectral curve can actually be realised as a solution 

to the discrete Bethe equations, a property which can be referred to as 
stability. In order to study the issue of stability, we find and explore the 
general two-cut solution or elliptic curve. It turns out that the moduli 
space of this elliptic curve shows a surprisingly rich structure. We 
present the various cases with illustrations and thus gain some insight 
into the features of multi-cut solutions. It appears that all admissible 
spectral curves are indeed stable if the branch cuts are positioned in a 
suitable, non-trivial fashion. 
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1 Introduction 

The Heisenberg magnet fl] is one of the very first quantum mechanical models, and it 
serves as the prototypical spin chain. Although it was set up almost 80 years ago it still 
remains a fascinating subject with many features left to be understood. Only three years 
after the discovery of the model and owing to its integrability, Bethe was able to write 
a set of equations 

f- =11 '- : for A; = 1,...,M, 

JT^k 
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which determine the complete spectrum |2|. The terms "Bethe equations" or "Bethe 
ansatz" later became synonymous for the exact solution for generic integrable spin chain 
models. 

Although the Bethe equations describe the complete and exact spectrum, it is virtu- 
ally impossible (and perhaps not very enlightening) to solve them concretely for generic 
states somewhere in the middle of the spectrum. Nevertheless some corners of the spec- 
trum are accessible (and interesting). In particular these are the low-energy and high- 
energy states for very long chains (the thermodynamic limit), where the Bethe equations 
are approximated by integral equations. Most studies have focused on the low-energy 
spectrum of the antiferromagnet and this regime is well-understood, see |3| for a review. 
For example, the antiferromagnetic state is a solution to an integral equation |4| and its 
excitation quanta are called spinous |5|. 

Conversely, the low-energy regime of the ferromagnet (or equivalently the high-energy 
regime of the antiferromagnet) is much less explored. The ferromagnetic ground state 
coincides with the vacuum of the Bethe ansatz, it is trivial. The excitation quanta are 
called magnons, and each magnon corresponds to a single Bethe root u describing the 
rapidity of the magnon. Magnons can form bound states which are usually called Bethe 
strings. In a fc-string centred at rapidity u there are k Bethe roots Mi,...,a; arranged in 
a regular pattern Uj = u + ^{k + 1 — 2j), i.e. the separation of adjacent constituent 
magnon rapidities is i, see Fig. ^, Of course this distribution pattern of Bethe roots 
is not exact when the length L of the chain is finite and in fact large deviations are 
observed. Nevertheless the string hypothesis can be used to perform a counting of all 
states which gives the expected exact result even at finite L \2\ (see also |6| for a recent 
account including references). 

An interesting type of "Bethe string" which has not found much attention until 
recently is one where the rapidity u scales as the length L of the chain. Here one 
distinguishes between short and long strings where the number of constituent magnons 
is either of 0{1) or of 0{L). For these strings the regular pattern is violated strongly, 
the distance between adjacent Bethe roots deviates much from i. Long strings were first 
investigated by Sutherland in [7J. This particular string consists of a condensate core 
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Figure 1: A Bethe string ' (a) [ and a short string (b) Bethe strings have a strictly 
regular pattern and only exist on infinitely long chains. For short strings, the 
separation of adjacent roots is of order 0{yV)\ the string shown here lives on a 
chain of length L = 10, 000. 
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Figure 2: Some long strings with mode number n 
a = MjL (L = length, M = number of roots). 
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and two tails, see Fig. 2b, Similar long strings were later considered in |8|. The distance 
of Bethe roots in the core is very close to i while in the curved tails the density of roots 
decreases to at the ends. Short strings were considered in |9, lO]. Unlike for standard 
strings the distance of Bethe roots is of 0[yL\ see Fig. lb, Short strings can as well 
be considered as very short versions of long strings |11|. The distinction of the two 
types is nevertheless useful because of their different energy scales in the limit of large 
L. Furthermore, long strings can be considered as smooth classical objects for which 
the discrete magnon constituents play a minor role. Short strings on the other hand are 
quantum objects and the number of constituent magnons must be a positive integer. In 
fact a short string is better viewed as a coherent state of bosonic magnons. In the strict 
thermodynamic limit the magnons do not interact with each other, their rapidities are 
merely influenced by the long strings which do interact non-trivially among themselves. 

Interest in the ferromagnetic regime was recently sparked by the AdS/CFT cor- 
respondence. Minahan and Zarembo showed the equivalence |10| of the spectrum of 
planar one-loop anomalous dimensions in the su(2) sector of A/" = 4 supersymmetric 
gauge theory to the spectrum of the Heisenberg magnet with zero overall momentum. 
Consequently it was observed that the energies of certain long ultra-relativistic string 
configurations in M x S"^ |12| agree precisely with the energies of certain long strings of 
the Heisenberg ferro magnet in the thermodynamic limit |8|, see |13| for reviews. This 
extended the earlier observation of |14| that the spectrum of magnons or short string 
excitations above the ferromagnetic vacuum agrees with the spectrum of quantum ex- 
citations of a short ultra-relativistic string orbiting the equator of S"^ . Some time later 
Kruczenski related both the ultra-relativistic limit of strings on M x 5^ and the ther- 
modynamic limit of the Heisenberg ferromagnet to one and the same Landau-Lifshitz 
model on S"^ |15|, see also |16|. 

The general solution for this classical model was constructed in |17, 18 1 in terms of 
spectral curves, see |19|. It was furthermore shown that the moduli space of spectral 
curves of any given genus has the expected dimension. The features of the spectral 
curves are in one-to-one correspondence to the integral equations obtained from the above 
thermodynamic limit of the Bethe equations. It is therefore clear that any state of the 
Heisenberg ferromagnet in the thermodynamic limit is approximated by a spectral curve. 
Less obvious is the question if every admissible spectral curve also has a corresponding 
solution to the Bethe equations, and if not, what are the stability criteria? In |20| one 
such criterion was derived from the requirement that the mode number of a long string 
can be defined unambiguously and self-consistently. The statement is that the directed 
density p of Bethe roots should not encircle the points in^ n ^ TL\ {0} when moving 
along the contour of a string. This is for example achieved if the density is bounded by 
\p{x)\ < 1 or |Imp(a;)| < 1, i.e. the distribution of Bethe roots should not be denser 
than the ideal Bethe string. In particular we expect \p{x)\ < 1 where the cut crosses the 
real axis. It is in fact a natural bound for the quantum model: The values p = in (along 
the imaginary axis) are distinguished because the pairwise interactions of Bethe roots in 



(1.1) become singular at that point. For the strictly classical model on the other hand 
the quantity p does not have a meaning, and solutions with density bounded by one are 
not at all different from solutions where the density exceeds one. A notion of stability 
exists in the semiclassical theory |21|; it demands the absence of tachyonic fluctuations. 



This condition does not coincide with \p\ < 1, it turns out that the latter is a stronger 
requirement in the case of a single cut |22|. 

So the question how to properly define stability remains. In this paper we would 
like to gain further insight into this issue by studying the general one- and two-cut 
solutions. The moduli space of two-cut solutions has two discrete parameters and two 
continuous ones, and it is sufficiently rich to investigate stability. In particular, we want 
to understand the role of condensate cores in the context of stability. In order to compare 
our analysis of the one- and two-cut spectral curves to actual solutions of the discrete 
Bethe equations (1.1), we further develop a numerical method for the construction of 
such solutions with large numbers of Bethe roots. Up to now, only few solutions that can 
be compared to corresponding spectral curves have been constructed |8|, and they only 
have a comparatively small number of Bethe roots. The numerical solutions with finite 
numbers of Bethe roots represent quantum states on a chain of finite length and hence 
can also be used to examine finite-size corrections to the thermodynamic limit. For one- 
cut configurations these corrections were studied analytically in |20,23,22|, some more 
configurations of roots were analyzed analytically and numerically in |24|. In this work 
we examine numerically the leading order finite-size effects explicitly for configurations 
beyond the one-cut case. 

This paper is organised as follows. We start in Sec. ^ with a review of the general 
solution in the thermodynamic limit by means of spectral curves. In order to illustrate 
our procedure on a simple example we will first reconsider and discuss the general one-cut 
solution in Sec. ^ The corresponding construction within the Landau-Lifshitz model is 
given in App. |Xj Sec. |4] contains the derivation of the general two-cut solution and its 
properties. In Sec. "Fwe will continue by applying the two-cut solution to study the issue 
of stability. To that end we have to determine the physical shape of the branch cuts in 
various regions of the parameter space. Finally in Sec. 7 we test our predictions against 
solutions of the Bethe equations with large but finite length, constructed numerically, in 
order to substantiate our claims. We conclude in Sec. 



2 Spectral Curves for the Heisenberg Ferromagnet 

We start by reviewing the spectral curves which describe solutions of the Bethe equations 
(LI) in the thermodynamic limit |18j. 

2.1 Baxter Equation 

To derive the properties of the spectral curves it seems convenient to consider the Baxter 
equation as advertised in |25 

T{u)Q{u) = {u + {fQ{u -i) + {u- ifQ{u + i). (2.1) 

This equation is fully equivalent to the Bethe equations in the following way: Let Q and 
T be polynomials of degree M and L, respectively. Then the set of solutions to the Bethe 
equation is equivalent to the set of Baxter equations. The Bethe roots Uk are simply the 
roots of the Baxter Q-function Q{u) = Y[k=ii'^ ~ '^k)- Furthermore the transfer matrix 



eigenvalue T{u) encodes the momentum eigenvalue P and the energy eigenvalue E in 
two equivalent ways as follows 

^" ^' - exp{+iP + iuE + C(m2)) + C»(u^), 



Tiu-i) 



iu — i)^ 



exp{-iP-iuE + 0{u^)) +0{u^). (2.2) 



In terms of the Q-function the momentum and energy read (unless L < 2) 



,P_Q(-i) ^_^Q'i+-2) ^Q'{—2. 



Q(+|) Q(+|) Q(-|) 

Note that for gauge theory and AdS/CFT the states are required to be cyclic, i.e. the 
net momentum must be zero, e^^ = 1. Here we will however not restrict to cyclic states 
but consider states with arbitrary momentum. 

Finally, the number of magnons M can be read off from the transfer matrix eigenvalue 
expanded around u = oo 

^ = 2 + (|L - i(L - 2M)(L - 2M + 2)) \ + Oil/u'). (2.4) 

Note that \{L — 2M)(L — 2M + 2) is the eigenvalue of the quadratic Casimir for a 
representation with spin ^L — Mr\ 

2.2 Thermodynamic Limit 

Before we take the thermodynamic limit, we make the following substitutions: Define 

exp(zp(x)) =exp(«/2x)-— — — and t{x) = , (2.5) 

where p{x) is called the quasi-momentum and t{x) is the properly rescaled transfer matrix 
eigenvalue. It is clear that the function p{x) has logarithmic singularities with opposite 
prefactors at x = {uk ± |)/L and a pole with residue i/2 at x = 0. Furthermore, we fix 
the ambiguity of p{x) by shifts of 27r by setting p{oo) = 0. Finally, t{x) is a polynomial 
in 1/x of degree L, i.e. it has an L-fold pole at x = and is analytic everywhere else. In 
these variables the Baxter equation reads 

t(x) = (l + i/2Lx)^exp ( , r ) exp(+ip(x - i/2L)) 

\ X — 1/2LJ 

+ (1 - i/2Lxfexp ( + ^^.^. ^ ) exp(-ip{x + i/2L)). (2.6) 

\ x + t/2LJ 

^Note the ambiguity M ^^ L + 1 — M which is related to the existence of mirror solutions of the 
Bethe/Baxter equations with M > ^L. These solutions have zero norm and thus are unphysical. Let 
us here restrict to physical states with M < \L. 



It is now easy to take the thermodynamic hmit L -^ oo. In this limit the magnon 
number is assumed to scale like M = La with fixed total filling a. The magnon rapidities 
scale like Uk ^ L and are assumed to distribute smoothly along certain contours in the 
complex plane with distance 0{L^). The distance between adjacent magnon rapidities 
defines the density p(x) of Bethe roots 

p{x) ~ for Uk ~ Uk+i ~ xL. (2.7) 

Mfe+l - Uk 

Note that the contours typically arrange vertically in the complex plane (but not neces- 
sarily strictly along the imaginary direction). Therefore the density function p(x) is in 
general complex and its phase is inversely related to the direction of the contour at the 
point X. For definiteness, let us decide that cuts generally go upwards in the complex 
plane, ImM^+i > Imw^. That means that the density will typically have a negative 
imaginary part, Imp(x) < 0. 

The Baxter equation in the thermodynamic limit becomes simply 

t(x) =2cos(p(x)). (2.8) 

The logarithmic singularities in p{x) move together to form linear discontinuities. These 
discontinuities are interpreted as branch cuts connecting various Riemann sheets oip{x). 
Apart from these and the pole at x = with residue i/2 the quasi-momentum p{x) 
is analytic everywhere. The function t{x) is also analytic except for an exponential 
singularity at x = originating from the pole of degree L ^ oo. Solutions p{x),t{x) to 
the equation (2.8) with the above properties describe the spectrum in the thermodynamic 
limit. Let us now study the properties of the solutions in more detail. 

2.3 Branch Cuts 

The function p{x) has discontinuities along certain branch cuts whose union shall be 
denoted by C. 

Density. The discontinuity of a cut is proportional to the density p{x) of Bethe roots 
(2.7[) in the vicinity of the point x 



Pi^) = 7r^{pix-^)-pix + ^)), xeC. (2.9) 

2m 

We shall take e to be an infinitesimally small positive number, i.e. p{x + e) denotes the 
limiting value of p at x G C towards the right of the cut. Note that the combination 
da; p{x) must be real and positive which determines the direction of physical branch cuts. 
The integrated density along a connected cut Ck will be denoted by the filling ak 



ttfc = / dx p{x). (2-10) 

'Ck 



Standard Cut. The function t{x) must remain analytic across a branch cut oi p{x). 
Equation (2.8) tells us that p{x) can change sign and shift by a multiple of 2tc without 
causing a discontinuity in t{x). 

If the sign changes across a branch cut we thus have 

p{x + e) + p{x — e) = 27rnfc for x G Ck, (2-11) 

where rij. is the (constant) mode number associated to the connected branch cut Ck- 
Using (2^ we can relate the density to the quasi-momentum 

p{x) = — (p{x - e) - -KUk) = — Uuk -p{x + e)). (2. 12) 

m Txi 

This type of cut can end in a square-root singularity x^ of p{x) where it takes the value 

p{x^) = link 

p{x^ + e) = nuk + 0{\^) for a;* G ^Cfc. (2.13) 

Note that this singularity of p{x) is indeed compatible with analyticity of t{x) at the 
singularity, cf. (2.8). At x=k the branch cut can be oriented in three different directions: 
Consider the radial coordinates x = x^: + re**^. The combination dxp{x) ~ e^^'^^'^y/rdr 
must be real and therefore the three possible orientations are separated by 240° (where 
the full rotation is by 720° due to the square root singularity). 

Condensate Cut. Conversely, if the sign does not change across a branch cut, the 
quasi-momentum must obey 

p{x - e) - p{x + e) = 27m';, for x e €[. (2.14) 

This type of cut would be required to end on a logarithmic singularity which is not 
compatible with analyticity of t{x) in pTsl). Nevertheless such a can exist if it ends on 
other cuts: Indeed, we can view a logarithmic cut (J2T4j) as the union of two parallel 
standard cuts (2 .11) . Therefore a logarithmic cut C can split up into two parallel cuts 
Cl (left) and Cr (right) at some point, see Fig. 3, Evidence for this splitting is provided 
by some numerical solutions to the Bethe equations for small L |9,8|. Compatibility of 



roots (2.9) to p{x) = —in'j^ and thus a logarithmic cut has constant integral density and 



(2.14) with (2.11) requires n' = tzl — "n-R. The integer n'^. determines the density of Bethe 
roots (2.9) to p{x) = —in'j^ and thus a logarithr 
extends strictly along the imaginary direction. 

2.4 Observables 



According to (2.3) the momentum P and energy E = EL appear in the expansion of the 



quasi-momentum at x = 



p{x) = — + P + Ex + Oix^). (2.15) 

2x 
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Figure 3: Two standard cuts Cl and Cpt with mode numbers ul and ur whose 
physical contours (as determined by the condition that p{x) dx be real) meet 
at a common point. The standard cut contours start and end at square-root 
singularities, which are indicated by crosses. Beyond their common point, the 
cut contours can be placed on top of each other - the result is a condensate 
cut as shown in the right picture. The density pc on the condensate cut can be 
determined by placing the individual contours infinitesimally close and parallel 
to each other, as shown in the left picture. It turns out that the density on the 
condensate cut is constant and purely imaginary: Pc{x) = iijiR — n^), see ( |2.14| ). 
Consequently, the contour of the condensate cut must be vertical. 



The relation between t{x) and the momentum and energy is slightly trickier: The limit 
of the expansion (2.2) gives two asymptotic expansions at the exponential singularity of 

t{x) at 



t{x) exp(-i/2x) = exp{+iP + ixE + 0(x^)) + 0{exp{-i/x)), 
t{x) exp{+i/2x) = exp(-iP - ixE + 0{x^)) + C(exp(+i/x)), 



(2.16) 



where the terms 0{u^) in (2.2) are interpreted as exponential singularities. Putting the 



two expansions together and comparing with (2.8) is in agreement with (2.15). 

The total filling a = M/L is easily obtained from the definition (2.5') of the quasi- 
momentum 



p{x) 



2a 



2x 



Oillx" 



(2.17) 



The expansion of t{x) around x = cxd in (2.4) yields a compatible expression, which 
however does not fix the above sign 



t{x) 



4x2 



Oiljx^ 



(2.18) 




Figure 4: Cycles of the spectral curve dp. Shown are the two sheets dpi, dp2 
of the Riemann surface of the spectral curve dp. The two sheets are connected 
by cuts Cj. To each cut Ci is associated an A-cycle Ai and a B-cycle Bi. The 
A-cycle encircles the cut contour, staying on the same sheet, while the B-cycle 
starts at re = CO on one sheet, passes through the cut and ends at x = oo on the 
other sheet. 



One can also derive two useful relations between the total and partial fillings 

Y^c^k, (2.19) 



a 



and the overall momentum and the mode numbers 

P = y^27rnfcafc- 



(2.20) 



2.5 Spectral Curve 

For the differential dp{x) the analytic structure is somewhat simpler: The shifts of p{x) 
by constants which arise when passing through a branch cut are not seen in dp{x). The 
differential dp{x) thus has merely two Riemann sheets with opposite signs and it is called 
the spectral curve. In particular, the condensate cuts drop out in the spectral curve. Only 
the standard branch cuts are seen in dp{x); they change the sheet or equivalently flip 
the sign. 

The above parameters of the solution can be read off from the spectral curve. First 
we need to introduce A-cycles and B-cycles for the cuts. A-cycles wind around the cuts 
while B-cycles extend from x = oo through a cut and back to x = oo but on the other 



Riemann sheet, see Fig. [4|pp "or a normal cut the A-cycle does not intersect with any 
branch cut of the quasi-momentum. The A-period of dp{x) is therefore zero 

dp{x) = 0. (2.21) 

The B-period yields the mode number Uk of the cut 

dp{x) = 2Tmk. (2.22) 

To understand this we split up the contour into the two parts before and after crossing 
the cut at point x. The first part yields p{x + e) — p{oo) = p{x + e) due to our normal- 
isation p{oo) = 0. On the other sheet the sign of the quasi-momentum is flipped and 
consequently the second part yields —p{oo) +p{x — e) = p{x — e). Altogether we obtain 
p{x + e) + p{x — e) and by (2.11) this equals 27rnfc. Note that a vanishing A-period of 
dp enables us to be rather unspecific about how the B-period returns to x = oo on the 
second sheet. 



Finally, the partial filling (2.10) of a cut can be expressed through the A-period of 

X dp{x) 

ak = I dx p{x) = i) dxp{x) = (h xdp{x). (2.23) 

This follows by substitution of (2.9) and partial integration. Note that the partial inte- 
gration assumes that the A-period of dp is zero. 

The above discussion shows that the cut contours and corresponding cycles can be 
deformed continuously without affecting the parameters of the solution. However, special 
care has to be taken when a cycle moves into a cut or passes through a singularity. In 
particular, the reality condition dp{x) > appears to play no important role in the 
purely classical approximation. 

In the case of many cuts the assignment of cycles is not necessarily unique anymore. 
Let us discuss what happens when two standard cuts join to form a condensate cut 
with four tails. Now there are various ways in which the four branch points could be 
connected, see Fig. 5, The standard method to set up the cuts is to connect complex 



conjugate pairs of branch points more or less directly, as in Fig. |5aj Then the above 
assignments of parameters works well. All cuts are of the standard kind and condensate 
cuts arise from two standard cuts running parallel for a while, see Fig. 3, 

Another option would be to connect the branch points horizontally with a condensate 
cut forming between the two standard cuts (Fig. |5b) . This option leads to non-zero A- 
periods 

(f dp{x) = ±2nn'-, (2.24) 



^Note that the B-cycles defined here are not cycles in the strict sense because they are open curves. 
Conventionally, B-cycles are defined as closed curves that pass through two branch cuts. A spectral 
curve with K cuts has genus K — 1 and hence only A' — 1 A-cycles and K — I closed B-cycles are 
independent. The set of K open B-cycles as defined here is equivalent to the set oi K — I closed B-cycles 
plus one open B-cycle. 

10 




(a) 



(b) 



(c) 



Figure 5: Possible cut contours (blue/red, tliick) and corresponding cycles for 
a core with four tails. In the standard configuration (a) , complex conjugate 
pairs of branch points are connected more or less directly. In this case, the 
A-periods vanish and the B-periods are non-ambiguous, even when the central 
cut segments join and form a condensate core as in Fig. 3. When the branch 
points are connected horizontally (b) , the mode numbers (B-periods) of the cuts 
are ambiguous, hence the configuration requires a condensate (black, vertical) 
between the two cuts. This implies that the A-periods do not vanish. Branch 
cuts that wind around each other (c) also require a condensate cut, hence complex 
conjugate branch points carry different mode numbers. 

where n'^ is the density of the condensate cut being intersected by the A-cycles. The non- 
trivial A-cycles then lead to ambiguities in the definition of mode numbers Uk and fillings 
ak- Also non-technically it is unclear how to associate suitable Uk and ak to these two 
individual cuts. The situation becomes even worse if the cuts wind around the branch 
points (Fig. 5c). Therefore the standard straight connection of complex conjugate pairs 
of branch points appears best to describe the parameters. 



2.6 Finite Gap Solutions 

The simplest types of spectral curves with the desired properties are the finite-gap solu- 
tions. These have finitely many branch cuts ( "gaps" ) and therefore finitely many branch 
points. The general ansatz for p'{x) is given by 

p [x) — 



x^h{x) ' 
where g{x) and h{x) are polynomials of degree K and 2K, respectively 



9[x) 



K 

E 



2K 



CkX 



K^) = W{X- Xk)- 



(2.25) 



(2.26) 



fe=i 
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The roots Xk are the square root branch points which must come in (complex conjugate) 
pairs. A pair is typically connected by a branch cut and therefore K is the number of 
standard branch cuts. The genus of the algebraic curve p'{x) equals g = K — 1. 

2.7 Stability 

Stability addresses the question which classical spectral curves with real and positive 
densities on the cuts can be realised approximately by solutions of the Bethe equations 
with large but finite L. An important stability criterion was derived in [20j 

^^^^n + tpix) ^ ^ forallnGZ. (2.27) 

Cfc n- ip{x) 

The criterion implies that the branch of the above logarithm can be uniquely defined on 
an isolated standard branch cut Ck- The argument of the logarithm 

n + ip{x) n{uk+i - Uk) + i Uk+n -Uk + i ,_ ^^. 

r^ ~ —, ^ : ~ : (2.28) 

n - ip{x) n[Uk+i -Uk) -I Uk+n - Uk - i 

approximates the scattering term in the Bethe equation (1.1). In a logarithmic form of 
the Bethe equations the condition implies that mode numbers can be unambiguously 
associated to the cuts, even in the full quantum theory. Unlike the classical quantities, 
here the density appears independently of the differential dx and therefore it is crucial 
to use the physical contour of the cut Ck- 

The most interesting case is n = 1; the conditions for n > 1 appear to be less 
constraining and perhaps redundant. It is quite clear that if the absolute density is 
bounded by unity everywhere, 

|p(x)|<l, (2.29) 

the stability criterion will be satisfied. However this condition is too restrictive in general. 
Nevertheless, for a stand-alone cut the density is typically highest at the centre Xq where 
the cut crosses the real axis and where the density becomes purely imaginary. It is then 
necessary to obey 

|p(xo)| < 1, (2.30) 

in order to satisfy the stability condition. Otherwise the argument of the logarithm in 
(2.27 ) will cross the negative real axis at Xq and wind around the origin once. 

What remains obscure at this point is how to interpret the stability condition for 
two standard cuts joined by a condensate. Also it is not clear if (2.27) alone can ensure 
that a given configuration of cuts can be realised by a concrete solution of the Bethe 
equations. In the following we shall study the general one-cut and two-cut solutions in 
order to shed some light on the various classically allowed configurations of cuts and 
when the stability condition is satisfied. 

3 One-Cut Solution 

In this section we review the simplest type of spectral curve with one cut. It has two 
parameters, the mode number n and the filling a. Its genus is zero and thus we will 
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only encounter algebraic and trigonometric functions. We also review the corresponding 
construction for the Landau-Lifshitz model in App. [Aj 

3.1 Solution 

The general one-cut solution was obtained in |18|. It has one mode number n, one filling 
a and it takes the form 



1 — ^-nnx \ Svmax 

p(x) = t™ H 4/I + - r- . (3.1) 

^^ ' 2x y (1 - 27rnx)2 ^ ' 

The derivative of the quasi-momentum defines the spectral curve 

p'{xf = (1 - 2^^^(1 - 2^))' (3.2) 

4a;^((l — 2nnxy + Snnax) 



and matches with (2.25[2.26[). The quasi-momentum is in agreement with the expansions 



at a; = 0, cxD (j2.15[2.17) and with the condition for the branch points (2.13). The higher 



terms in the expansion at a; = yield the total momentum and energy 

P = 27ma, E = 47rVa(l-a). (3.3) 

If we choose to restrict to cyclic states as required for AdS/CFT we have to set P = 27rm 
with integer m. Then the total filling must be a rational number a = m/n. However, 
all these solutions are unstable. As we shall see later this is related to the fact that the 
total momentum leaves the first Brillouin zone, \P\ > vr. 

3.2 Cut Contour 

The physical contour of the cut is not a simple function. In particular, it is not given 
by the natural branch cut of the square root in p{x); it lies somewhat closer to the 
origin. In order to find the contour, we shall make use of the identity (2.12) which 
relates the density to the quasi-momentum. The integrated density must be real and 
positive; therefore we will need the integral of the quasi-momentum 



, 1 — 2TTnx / Simax 

A{x) = / p[x) ax = nnx H w 1 



;i - 27mx)2 

1 — 2Tmx I Siinax 

artanh 4 / 1 + 



1 + 2TTnx y (1 — 2-Knxy 



, l — 2Txnx / Svrnaa; , ,„ ^, 

aarcoth ^ / 1 + r- . (3.4) 

' 1 - 2Tinx - 2a y (1 - 2TinxY ' ^ ' 

On the cut the combination A{x) — nnx must be real. At the branch points 

l-2a±2V«(l-a) 
"^ = 2^ ('•') 
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Figure 6: Physical position of the cut (sohd/red) for a one-cut solution with n = 
2, a = 0.05. Two alternative contours with positive filling originating from the 
same branch points are shown: The second (dashed/blue) contour corresponds to 
the conjugate one-cut solution with n = —2, a = 0.95. This contour is unphysical 
because the filling exceeds 0.5. The third contour (dotted/green) extends to 
infinity, it is inconsistent because the filling is infinite. The remaining lines depict 
possible cuts with real filling not ending on the branch points. Fluctuation points 
are marked by dots. 



it takes the values j^(x±) 



Tmx± = ±|7ra and thus the contour is defined by A{x) —imx G 
^Tra[—1, 1]. In Fig. "6" we have plotted the contour of a sample cut. 

In fact as discussed below (2.13), a branch cut with positive density can, in principle, 
originate from a branch point in three different directions. The shortest path will turn 
out to be the only physical choice. 

The longer circular path which encircles the origin is equivalent to the one-cut solution 
with opposite mode number n' = —n and conjugate filling a' = 1 —a. This can be seen as 
follows: We deform the shorter cut continuously to the longer cut by rotating it towards 
the right by 240°.^ At some point the cut has to pass the point a; = oo. This has the 
following two effects. Firstly, the A-cycle describing the filling (2.23) intersects a; = oo 
twice. 



This adds to a the residues | 



circulation) so we obtain a' = a + 2{ 



a and ^ 



a (the opposite sheet with opposite 
a) = 1 — a. Secondly, the point x = oo now 
resides on a different sheet. This implies p(oo) = 27rn and we have to subtract the 
constant 2?™ to normalise to p(oo) = 0. At the branch point that leads to p{x^) = — vrn, 
i.e. the new mode number is n' = —n. Note that this configuration has a large filling 
a > ^ because the shorter cut has a < |. It is therefore unphysical. 

For the third choice we rotate the contour towards the left by 240°. This contour 
escapes towards ±ioo on both sides. Moreover the density p on the contour approaches 
a finite value at infinity. Therefore the total filling on the contour is infinite, and the 
configuration is inconsistent. 

^Alternatively one can rotate by 120° and flip the sign of p. 
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Figure 7: Absolute density |/9(xo)| at tlie centre of a single cut with mode 
number |n| = 1 versus the cut's filling fraction a (blue/thick). Also displayed are 
the lines 1/m for integer m. According to equation (3.7), the intersection points 
of these lines with the density mark the maximal fillings acond of stable solutions 
with mode numbers m (see also Tab. IJ. 

3.3 Stability 

We are now ready to consider the issue of stability. The stability condition for the one- 
cut solutions implies that the density must be bounded by unity where the cut crosses 
the real line. Assume the cut crosses the real axis at xq which must be the solution of 
the equation A{xq) = nnxo. This equation is transcendental and we can only solve it 
numerically for given values of n and a. Once we have the value xq the absolute density 
is given by 

\p{xo)\ = -\p{xo) -7in\. (3.6) 

n 

The density for solutions with n = 1 is plotted in Fig. 7. We are interested in the filling 
ttcond where the density reaches the maximum allowed value, |p(a;o)| = 1 or p{xo) = 
7r(n + 1). Let us assume for simplicity that n > 0. Because of the way p{x) and A{x) 
scale with x and n, the intersection point of a cut with mode number n and filling a and 
the density at that point are directly related to the intersection point and the density of 
the cut with mode number 1 and the same filling: 



Xo{n, a) 



xo(l,a) 



n 



pn{xo{n, a)) = npi{xo{l, a)) . 



(3.7) 



Hence it is sufficient to obtain the fillings a at which the absolute density at the centre of 
the cut with mode number 1 equals 1/n. These fillings are exactly the maximal fillings 
of the cuts with mode number n that are allowed by the stability criterion; they are 
indicated in Fig. |7| 
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n\ 


Q^cond 


^cond/47r2 




\n 


C^cond 


^cond/47r2 


1 


0.2092896452 


0.1654874896 


6 


0.0069071371 


0.2469394280 


2 


0.0596024470 


0.2241999811 




7 


0.0050818745 


0.2477464054 


3 


0.0271903146 


0.2380590127 




8 


0.0038944180 


0.2482720935 


4 


0.0154373897 


0.2431852264 




9 


0.0030790284 


0.2486333844 


5 


0.0099228217 


0.2456089819 




10 


0.0024951483 


0.2488922545 



Table 1: Maximal fillings acond and the corresponding energies -Econd for one- 
cut solutions for the first few mode numbers n. Above these fillings, solutions 
require a condensate for stability. 



Alternatively, one can solve the equation p{xq) = 7r(n + 1) for xq in terms of a and 
n. The solution is substituted in A{xq) = imxo which yields the equation 

•^^ ^1 — g^ — 2nq'' 



2{q + n) arcoth(g + 2n) 



1-q' 
2n 



artanh 



2n 



(3.8) 



Here q is an auxiliary variable that encodes the maximal filling ctcond and the intersection 
point xq 



q^ 



1 



ttcond = -r~i — \ — V ' ^0 = 7w — \ — V ■ (3.9) 

4n(g + n) 2TT[q + n) 

We list the first few values of acond in Tab. 1, Note that ctcond < l/2n, which implies 
that stable solutions have a momentum P < n. Hence, cyclic one-cut solutions cannot 
be stable. In fact one can solve the above equation perturbatively using that q ~ 1/n. 
The expansion of Ocond reads 

1 7 367 540373 1895953 ,^, , ^o, , 

+ 0{l/n^^). (3.10) 



C^cond 



+ 



+ 



4n2 144n4 38880n6 293932800n8 5290790400n 
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3.4 Fluctuations 

Fluctuations are very small cuts which can exist in a background the one long cut |8|. 
For the one-cut solution they were discussed in |8,26|,|22j. Their position x^, is determined 
by the long cut, 

p{x^,) = vrn*, (3.11) 

but they are not strong enough to back-react on the position of the long cut substantially. 
We shall again assume that n > 0. The solution to the above equation reads 



— = 2?™ 1 - 2a) + 2TT{n, -n)Jl -. — ^ —^ 

X* y (n* - n)-^ 



(3.12) 



The momentum and energy of a fiuctuation mode quantum with 6a = 1/L are given 
by [26^ 22] 



6P 



2?™* 



6E 



Att' 



n{2n^ - n)(l - 2a) + {n^ - nf 



An^ail — a) 



(n* -ny 



(3.13) 
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If the long cut is rather short, the fluctuations will reside close to their vacuum 
positions x* ~ l/27rn=K. For fluctuations of the same mode number as the cut, n^, = n, 
the condition p(x*) = vrn,,, is solved by the branch points, which means that this particular 
fluctuation will increase the size of the long cut. The longer the long cut gets the more 
will it attract fluctuations with nearby mode numbers n^ ^ n towards it |22|. At some 
value of a, the fluctuation with n,, = n + 1 collides with the long cut (from the left). 
This happens precisely when the density reaches unity at the real axis, cf. (|3.6|, and so 
the above discussion applies. Note that at this point the fluctuation with n* = n — 1 is 
still at some distance from the cut. 

In conclusion, we can infer that for all stable one-cut solutions the fluctuations are 
well-separated from the long cut. When the fluctuation collides with the cut we can 
actually argue independently for an instability: Now the long cut can be filled not only 
from both ends, but also from the middle. In practice we expect that a condensate cut 
(with unit density) will form in the middle of the long cut. 

3.5 Condensate Formation 

Let us add a vertical condensate cut ending on the existing branch cut, see Fig. 8, This 
is achieved by shifting the quasi-momentum p{x) by — 27r in the region enclosed by the 
two contours, cf. Fig. 9, Clearly the new curve satisfies all conditions for the classical 
spectral curve. In effect, the shift decreases the mode number on the inner part of the 



original branch cut by one unit according to (|2.11). Furthermore, the density function 



is changed according to (2.9,2.12) and thus the inner part of the branch cut has to be 
moved to obtain a real density 

p(x) = — (7r(n-l) -p(x + e)). (3.14) 

TXl 

The condensate cut can end on any point of the branch cut, and we show the various 
potential configurations in Fig. 8, Any of these configurations with a condensate appear 
to be possible from the point of view of spectral curves. However, the Bethe equations 
will single out one particular configuration as the distribution of Bethe roots in the 
thermodynamic limit. To understand the physical distribution we have to distinguish 
three qualitatively different cases, see Fig. 8, 

In the first case the filling is below the threshold for condensate formation, a < 
ttcond, as discussed in the previous subsections. Here all potential deformations leave the 
branch cut almost vertically and then circle around the origin. Due to the residue at the 
origin the filling of the configuration is altered and becomes larger than ^. Inserting a 
condensate cut therefore leads to an unphysical configuration when a < Ocond- 

At a = Ocond the fluctuation point with mode number n* = n + 1 crosses the branch 
cut and effectively acquires the mode number n^, = n — 1. We now have two fluctuation 
points with the same mode number n — 1. For definiteness, let us call the new point 
n^, = {n — 1)'. When increasing a further, it will eventually collide with the other 
fluctuation point. This happens at a = Ocrit with 



11/1 , , 
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(a) a < a. 



cond 



(b) Qfcond < a < Qfcrit 



(c) acrit < a 



Figure 8: Potential deformations (thin/black) of a brancli cut (solid/red) after 
inserting a condensate (dashed/green) for n = 2 and[(a)JQ = 0.03, (b) a = 0.065, 
'(c)' a = 0.8. The physical deformation (solid/blue) exists only for a > Ocrit- The 
contours originating from near the ends of the branch cut close up after encircling 
the origin. 






Figure 9: Possible deformation of a branch cut: For large enough filling a > 
Ocondi a closed loop cut (thick/green) can be added to the branch cut (thin/red). 
In the region enclosed by the loop cut, the quasi-momentum p{x) is changed to 
—p{x) + 27rnL, where til is the mode number of the loop cut. As shown in Fig. 3, 
the density pc on the common part of the two contours equals i{ni^ — ti-b), where 
nB is the mode number of the branch cut. Hence the resulting configuration has 
a condensate core, as shown in the right Figure. 
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(a) 



(b) 



(c) 



(d) 



(e) 



Figure 10: For large fillings a, one-cut solutions are degenerate cases of two- 
and three-cut solutions. When a macroscopic number of Bethe roots gets added 
to the fluctuation point (blue dot) in (a) , two square-root branch points form, 
the closed loop cut splits at its cusp and yields a genuine two-cut solution (b) , 
This singles out the deformation shown in (a) as the physical one for the case 
Ocond < a < Ocrit- At even higher fillings a > acrit (c) an excitation of the 
fluctuation points results in a three-cut solution (d). A local minimum of the 
energy is reached when all roots from the third cut are moved to the second, 
leaving a two-cut configuration and a bare fluctuation point as in (e) , 

Consider now the case acond < a < Ocrit- The deformations originating from the 
branch cut can now have two qualitatively different shapes. The contours originating 
from near the ends of the branch cut behave like in the first case above and are thus 
unphysical. The contours originating from near the centre however are just small defor- 
mations of the original cut and consequently they lead to the same filling. All of these 
configurations are in principle okay, however, only one can be physical. Indeed one of the 
configurations is special: the limiting shape, i.e. the largest possible small deformation. 
It is distinguished by a cusp at the fluctuation point for n,,, = (n — 1)', see Fig. 8b, We 
can argue that this is the physical configuration: When we add macroscopically many 
Bethe roots to the fiuctuation point, it will split up and form two square root singu- 
larities. The deformed cut will split at the cusp and form a condensate with four tails 
as in Fig. |10b[ This configuration is a genuine two-cut solution as we shall see in the 
next section. When we take the new Bethe roots away we should return to the one-cut 
solution. This is possible only if the deformed contour meets the fiuctuation point with 
mode number n^ = [n — 1)', as in Fig. 8b 10a, In conclusion the one-cut solution for 
ttcond < a < ctcrit IS a degenerate case of a two-cut solution. 

The final case is Ocrit < «• Here the two fiuctuation points with mode numbers 
n=K = (n — 1)' and n^, = (n — l) have joined and branched off into the complex plane. The 
configurations are similar to the above case of Ocond < a < Ocrit, but now the limiting 
shape meets both fiuctuation points and thus has two cusps joined by an approximately 
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vertical contour. In this case we cannot add Bethe roots to any of the two fluctuation 
points individually because it would violate the reality condition (reflection symmetry of 
the conflguration about the real axis). They can only be excited in pairs to create four 
new square root singularities and a three-cut solution as in Fig. lOd, In other words, the 
conflguration is a degenerate case of a three-cut solution. The third cut is the line segment 
joining the two fluctuation points. By removing Bethe roots from this cut one lowers the 
energy and we might therefore call this solution unstable. Nevertheless we expect it to 
be a perfectly well-behaved solution of the Bethe equations in the thermodynamic limit 
albeit with non-minimal energy. A local minimum is obtained by shifting all Bethe roots 
from the third cut to the second, as in Fig. lOe. This is a true two-cut solution which 
will be discussed in generality in the next section. 

Finally we remark that for n = 1 the third case does not exist because Ocrit = |- 
This is because the fluctuation point with n,,, = is always at x* = oo and cannot join 
with the one for n^ = 0'. Therefore the solution with n = 1 represents a local minimum 
of the energy for all < a < 1/2. 

4 The General Two-Cut Solution 



As described in Sec. 2, a solution to the Bethe equations (1.1) in the thermodynamic 
limit is given by a quasi-momentum p, whose domain is a multi-sheeted cover of C. The 
corresponding spectral curve dp has B-periods in 27rZ and vanishing A-periods, while 
the quasi-momentum p{x) has a simple pole of residue 1/2 at a: = 0. In this section, 
the general spectral curve with two cuts will be constructed and its properties will be 
investigated. Spectral curves with two cuts have a Riemann surface of genus one and 
thus can be described in terms of complete elliptic integrals of the flrst, second and third 
kind. These are denoted by K, E and 11 and are deflned as 

K(g) := / ^dt, 

Jo v/(l-t2)(l-gt2) 





^<''-' ^^*- 



n{z\q):=[ ^dt. (4.1) 

The two-cut solution can be obtained by direct integration of the general ansatz 



(2.25). Here a different approach will be followed which makes use of the result obtained 



for the symmetric case in |27|. 



4.1 The Symmetric Two-Cut Solution 



We will obtain the general two-cut solution by generalising the symmetric solution [28,27 

poW - -^^J#l^n(^ ,). (4,2: 
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where q = 1 — a^/h"^- As a function of z, the eUiptic integral of the third kind T{{z\ q) 
has branch points at 2 = 1 and z = 00. Its discontinuity across the branch cut between 
these two points (from the lower to the upper half plane) equals 27ri times the residue of 
the integrand in (4.1) at t = 1/v^: 

Iliz + ie\ q) — Iliz — ie\ q) = (h dt 



■-^'v^W^y f^'^) 



The function n(z) := n(gz^/(z^ — a^)\ g) used in (4.2) is chosen such that it has branch 
points {ao, 60) ~c^0) ~^o}) hence T^{z) has two branch cuts that lie symmetrically around 
2 = 0. In principle, Oq and 60 can be arbitrary complex numbers. However, physical 
branch cuts must be symmetric about the real axis, therefore only the case 60 = clq is 
physical. The discontinuity across the cuts can be inferred from (4.3) and equals 




Ii{z + e)-Yi{z-e) = ±-. ";, " ' . (4.4) 



In (4.2), the factor in front of 11(2;) makes the discontinuity constant, but also switches 
sign across the cut. As a result, the sum of the limiting values of po on either side of a 
cut is 

Pq{z + t) +po{z - t) = ±'n/\n. (4.5) 

This means that for An G 2Z, the Bethe equations (1.1) are satisfied and the cuts (aoi ^0) 
and (— flo; ~^o) have mode numbers An/2 and —An/2 respectively.^ Since 11(0) = K(g), 
Po{z) has a simple pole at z = with residue 

Po(2) = -AnaoK(g)- + ... (4.6) 

z 

For hQ = (Iq, q = 1 — exp(4i argag) lies on a unit circle centred at 1 ^ and hence depends 
only on the argument of oq, not on its modulus. As a result, one real degree of freedom 
remains after imposing the condition resoPo(-2) = 1/2, which corresponds to the filling 
of the two symmetric cuts. 

4.2 Construction of the General Two-Cut Solution 

The solution (4.2) can be generalised to the non-symmetric case by forming the compo- 
sition pq o /i, where /x is a Mobius transformation 

u:C^C, xi— *-z = jiix) = . (4.7) 

rx + s 



*The sign in (4.4) depends on the choice of ag, 60 and e. For Re(ao) > 0, Im(ao) sufficiently small 
and e > 0, the positive sign holds on the left, the negative sign on the right cut. 

^Conventionally, the elliptic integrals have a branch cut on the real axis for q > I. The unit circle 
crosses the cut aX q — 2 and for Im q < one has to analytically continue the elliptic integrals suitably. 
Alternatively one can apply the Landen transformation used in |11 which maps the unit circle to the 
interval [0, 1] where no branch cut is encountered. 
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It turns out that transforming the symmetric solution this way preserves enough of its 
structure while providing sufficient freedom for constructing arbitrary two-cut solutions. 
By suitably choosing the Mobius transformation /z and the modulus g, one can map 
any four points {a,b,c,d} G C onto the branch points {ao,bQ,—ao,—bo}, i.e. one can 
construct a function with any four branch points {a, b, c, d} by solving the system of 
equations 

fi{a) = ao , fi{b) = bo , fi{c) = -oq , fi{d) = -bo (4.8) 

for ao, bo and the parameters of the transformation /i. In fact, the equations (4.8) do not 
completely fix oq, 6o and /i. However, only solutions whose cuts are symmetric under 
refiection about the real axis are physical. Hence it is useful to set 60 = '^0 and restrict 
oneself to real Mobius transformations (i.e. s, t, m, r e M), which map complex conjugate 
pairs to complex conjugate pairs. 

Applying the transformation to po directly would move the pole from z = to 
X = yU^^(O) = u/t. In order to have the pole at x = in the transformed function, it 
has to be moved to z = /i(0) = u/s in the original function pq. This can be achieved by 
adding a term to po 



al{bl~z^) al-z^ AnaoK{q)u 
Po{z) .- Poiz) - y ^2(^2 _ ,2) -^f- ^(^, _ ^) • (4.9) 

The prefactor is necessary for retaining the structure of cuts of the original function 
without introducing additional poles at z = ±ao. The function to be transformed now 
reads 



Po[z) 



An I alibi - z^) 



aoz{zs — u)\i bl{al — z"^) 
L{al-z^)K{q) + al{zs-u)Ii{^-^^\q)\ , (4.10) 

and the candidate for the generalised solution is 

p:=Poo^i + C. (4.11) 

The constant C must be chosen such that p{oo) = 0, therefore it has to be defined as 

C := -Po o /i(oo) = -Po it/r) . (4.12) 

The function p has branch points {a, 6, c, d} = /i~^({ao,6o, — c^o, ^^0}) and branch cuts 
Ci := (a, 6), C2 := (c, d). As can be verified directly by computing the derivative of (4.11), 
dp is of the form (2.25). The following sections will show that the function p indeed meets 
all requirements for being a valid quasi-momentum while providing maximal freedom in 
the choice of physical parameters. 
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4.3 A-periods and Mode Numbers 

For being a valid quasi-inomentum, p must have vanishing A-periods and integral mode 
numbers (B-periods), as was established in Sec. 2.5, Because the function (4.11) is 
single-valued by construction, the integrals of dp over A-cycles vanish, as long as one 
chooses the cuts Ci and C2 of p between the branch points (a, 6) = /i~^({ao, 60}) and 
(c, d) = /i~^({— oo, —bo}) in a way that they do not cross each other. Since p vanishes at 
infinity, its mode numbers ni and n2 are (cf. ( 2.22, )) 

27ml := Poix + e) + Po{x - e) 

= poifJ'ix + e)) + po{fi{x - e)) + 2C 

= 2C + 7rAn, x e Ci, (4.13) 

and similarly 

27m2:=Po{x + e) +po{x - e) = 2C -nAn, x e €2- (4-14) 

For obtaining physical solutions, the parameters of the function p must be chosen such 
that Hi and n2 are integer numbers. 

4.4 Energy and Momentum 

The residue of the pole as well as the total energy and the total momentum of a solution 
p are encoded in the expansion of p{x) at x = 0. Performing this expansion, comparing 
with (2.15) and demanding that the residue of the pole be 1/2 yields the equations 

u^ iplrs - tu) + als {ru^ + hls^^st - 2ru)) -,^, A „ 
^^ HQ) - ^^ 



E 



aou {st — ru) {h'^s^ — u 
An{st — ru)s'^W{u/s) 



■u^)(aQS^ — u^) 



hi E(g) - '^"""^ o7.L .T.^r^r ' K(?) I , (4.16) 
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4.5 Fillings 

The partial filling fractions ai, 02 can be obtained by finding the total filling a with the 
help of (2.17) , 

1 An f (airs — tu)(blrs — tu) ^^, , , , ,^, ,\ ,..^, 
a= ^^^^^ , ^ ^^ — , , ^^ ° '- Kg - hoist - ru) Eg , 4.18 

2 r^W(t/r)\ bo{st-ru) ^^' "^ ' ^^7' ^ ' 



and using a = ai + ^2 together with riiai + ^20:2 = P/2'k (2.20). 

Alternatively, ai and a2 can be calculated directly by performing the contour integra- 
tion (2.23) over p{x) along A-cycles around the cuts. This is done explicitly in App. C, 
the result is 

/\n ( hlnist-rul , , vr.^ (^^^, JoJ ' BiPl ' Vs) \ ,, , , 
"^'^ = h-^ [ 2rW(t/r) ^(^) + 2R37^ T ^"^'^ WF) ) ^'^ 

± f E(g) K(g) ^ ^ K^(g) ± fM^^ K(g) fl(./.) 
t ut OoU[st — ru) 

± b!:i^tll^ E(g) flit/r) ± -lr{alrs-tuMrs-tu) 

(4.19) 



Provided that resop(x) = 1/2 holds (which can be implemented by imposing (4.15)), one 
can verify that indeed ai + a2 = a and rtiai + 71,202 = P/'^'^i as required by (2.19) and 
(2.20). 

4.6 Solving for Physical Parameters 



The parameters of the solution ( |4.11 [) are oq, &0; An and the three parameters of the 
Mobius transformation yU. The solution p is only physical if these parameters are adjusted 
such that all physicality conditions are satisfied. Namely, the mode numbers (4.13) and 
(4.14) must be integral, the filling fractions (4.19) must be real and (4.15) must hold. 
The correct behaviour at infinity (2.17) is provided by construction. 

Setting 60 = (^0 and restricting to real Mobius transformations results in complex 
conjugate branch points, which guarantees that the filling fractions are real. Further, 
equation ( |4.15| ) can be solved for t (or r, alternatively): 

t= — - ^^ W{u/s) K(g) . (4.20) 

S Oq 

Imposing these constraints leaves the complex oq, the real An and two real parameters of 
/i as free parameters. Taking into account that the curve p is invariant under the rescaling 
{ao,t,u} -^ \{ao,t,u}, A G M, only four real parameters remain. They correspond to 
the freedom in the choice of the mode numbers ni, n2 and the filling fractions ai, 02- 
This is consistent with the result of |18j, according to which the algebraic curve with k 
cuts has 2k parameters: k discrete mode numbers and k continuous fillings. 
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The expressions for the mode numbers (4.13) and (4.14) can be easily solved for An, 
yielding 

An = ni - n2 . (4-21) 

The physicality constraints cannot be solved analytically for the remaining parameters, 
hence one has to resort to numerical methods for finding explicit solutions p for given 
mode numbers and filling fractions (or other sets of physical quantities, such as energy 
E and total filling a). 

4.7 Finding the Cut Contours and the Density 

For a given solution p, the physical contours of the two branch cuts are determined by 
the condition that p{x)dx be real. Together with the relation (2.12) between the density 
and the quasi-momentum, this condition yields a first-order differential equation for the 
cut contours. The physical cut contours can be obtained by numerically integrating this 
equation. Alternatively, one can obtain an explicit expression for the integrated density, 
similar to (3.4) in the one-cut case. The integral of the density reads 

' I r 1 

p{y)dy = — {p{y) -7mj)dy= —{A{x) - A{x^)) +inj{x-x^), (4.22) 

where A{x) is the integral of the two-cut quasi-momentum p{x) and rij is the mode 
number of the respective cut. The integral A can be calculated analytically and is given 
in App. D, On the physical contour, the expression (4.22) must be real when x=k is one 
of the branch points. 

Note that the function A{x) can also be used for obtaining initial positions of Bethe 
roots for the computation of discrete Bethe root distributions. Namely, when one wants 
to approximate a given two-cut spectral curve p with a solution to the discrete Bethe 
equations on a chain of length L, then this solution has to consist of M = aL roots. 
Mj = ajL of these roots have mode number nj and must lie on or close to cut Cj, j = 1, 2. 
Taking the solutions Xk to the equations 

1 A; — 1/2 
— {A{xk)-A{x,j))+inj{xk-x,j) = -^ , k = l,...,Mj, 3 = 1,2 (4.23) 

as initial positions for the Bethe roots reproduces the density function on the respective 
classical cut well and thus gives a good chance of numerically finding an approximating 
solution of the Bethe equations. Of course, the same method can be employed for the 
one-cut solution, using the respective integral (3.4) . This way, the solutions shown in 
Fig. lb, 2 were obtained. 



In Sec. 7.1, we describe another simple way for obtaining discrete distributions of 
Bethe roots which become exact in the limit of large length L and fixed number of roots 
M. They can be easily constructed for any number of cuts and we use them to obtain 
numerical solutions of the Bethe equations in Sec. 7.2 , However (4.23) gives a better 
approximation when L is not sufficiently large. 
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5 Moduli Space for Consecutive Mode Numbers and 
Stability 

In the last section, the general two-cut spectral curve has been constructed. For given 
mode numbers ni and 712, the moduli space of this solution is two-dimensional. It can be 
parametrised by the two fillings ai and 0:2, by the total filling a and the energy E or by 
any other pair of physical quantities. In this section, the moduli space of configurations 
with two consecutive mode numbers is investigated. It turns out that this space has a 
very rich structure and even is globally connected, that is the spaces of different pairs of 
mode numbers are interconnected. Some features of the moduli space for non-consecutive 
mode numbers are presented in the following section. 

5.1 General Features of the Moduli Space 

As discussed in Sec. 3, a single cut with mode number ric and a small filling is centred 
around l/(27mc). When the filling of the cut increases, the neighbouring excitation 
points with mode numbers nc± 1 get attracted by the cut, until the excitation ric + l (for 
positive ric) collides with the cut when the absolute density at its centre reaches unity. 
The same happens when the excitation n^ + 1 has a small but finite filling, i.e. when it 
is replaced by a small cut. As will be shown in more detail below, larger cuts in general 
attract smaller cuts when their filling is increased. 

To begin with, the moduli spaces of configurations with two consecutive mode num- 
bers can be divided into several regions, as is shown in Fig. 11||12[13[14|15 for the cases of 



rii = 1, n2 = 2 and rii = 2, n2 = 3. Solutions in each region and on each line separating 
the regions share certain characteristics. In the following sections, the different regions, 
lines and points of the moduli space will be discussed in detail and examples will be 
given. For simplicity, only non-negative mode numbers rij > are considered, which 
correspond to cuts that lie to the right of the imaginary line. The discussion for negative 
mode numbers is completely analogous. The two cuts are named "cut one" and "cut 
two". Unless otherwise stated, the mode number rii of cut one is smaller than the one 
of cut two, n2 = ni -t- 1. 

In all figures that show cut contours or densities, cut one has blue colour (thick lines), 
while cut two has red colour (thin lines). Plots of cut contours in the complex plane 
always have an aspect ratio of 1 : 1. In all density plots, the density is shown versus the 
distance from the cut's centre, measured along the cut's contour. 

5.2 Regions 

The moduli space of configurations with two given consecutive mode numbers can be 
divided into three regions. Configurations in each region share certain characteristics: 
In region I, the two cuts are disjoint; in region II, the cuts join with each other and form 
a condensate with four tails; in region III, the two cuts are disjoint but require a third, 
closed cut and a condensate for stability. In the following paragraphs, the characteristics 
of the different regions are described in more detail and example solutions are shown. 
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Figure 11: Moduli space of configurations with mode numbers ni = 1, n2 = 2, 
shown in the (a, E) plane. Regions of the space are named by roman numbers 
I-III, while lines are assigned Latin letters A-J. Special points have labels X- 
Z. The dotted area is magnified in Fig. 12. The moduli spaces of higher mode 
numbers look very similar, cf. Fig. 14 
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Figure 12: Magnification of the dotted area in Fig. 11 
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Figure 13: The same moduli space as in Fig. [ill shown in the (01,02) plane. 
The lines Ai, Bi, A2 and H2 lie on the axes. 
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Figure 14: Moduli space of configurations with mode numbers ni = 2, n2 = 3, 
shown in the (a, E) plane. The dotted area is magnified in Fig. 15. Structurally, 
the moduli spaces of higher pairs of mode numbers like this one are very similar 
to the case ni = 1, 77-2 = 2 shown in Fig. 11, One can already see that the part 



of the space shown here is connected to the part shown in Fig. 11 
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Region In,n+i'- Two Disjoint Cuts. In this region, the contours of the two cuts with 
mode numbers n and n + 1 are disjoint, and the cut with mode number n + 1 lies further 
left (closer to the origin) than the one with mode number n. The absolute density is 
bounded above by 1 everywhere on both cuts. An example is shown in Fig. 16, When 
the filling of one of the cut increases, this cut attracts the other one, as can be seen in 
Fig. pTpS 



Region Iln^n+i- Condensate with Four Tails. In this region, the standard con- 
tours^ of the two cuts bend towards and cross each other. An example is shown in Fig. 19| 
As explained below ( |2.14| ), the sum of the densities at the two common points of the 
standard contours equals —i (since the difference of the mode numbers is ni — n2 = —1) 
and solutions in this region allow for straight condensate cuts with constant density 
Pc{x) = —i between the two common points of the cuts. The condensate-cut version of 
the configuration in Fig. 19 is shown in Fig. 20. Beyond the two common points, the 
remaining parts of the standard cut contours form four tails that are attached to the 
ends of the condensate. The tails of cut two (with mode number n + 1) lie further left 
than the ones of cut one (with mode number n). 

Whenever a particular solution p allows for a condensate cut, it will be assumed that 
the configuration with the condensate cut is the physical version of the solution. This 
assumption may not seem very motivated at the moment, but will become more justified 
during the subsequent discussion. A first sign of its verity is the fact that there is no 
well-defined direction on standard cut contours that cross each other. This can be seen 



as follows: According to ( 2.12| , the density p on a standard cut contour is proportional 



to the difference between the limiting values of p{x) on either side of the cut contour: 
p{x) ~ p{x + e) — p{x — e). When a standard cut contour crosses another standard cut 
contour, the value of p{x ± e) on either side of the contour changes to —p{x ± e) -|- 27in, 
where n is the mode number of the cut that is traversed; this follows from the Bethe 



equations (1.1). Hence the density p on the contour switches its sign, and in order 
that the corresponding root density p{x) dx stays positive, the direction on the contour 
must change as well. However a unique direction on each contour was assumed in the 
derivation of the spectral curve and quasi-momentum in Sec. 2.3. This argument makes 
it seem very unlikely that configurations with crossing standard cut contours can be 
realised as limits of Bethe root distributions. 

As indicated by the dashed lines in Fig. |11|12|13[14|15| region II„,„+i can be further 
subdivided. Between the lines Dn,n+i, Gn,n+i and C„, the standard cut contours have 
the usual form, as in Fig. 19: Neither do they encircle the origin, nor do they extend to 
infinity. In contrast, beyond line Gn,n+i the standard contour of cut one encircles the 
origin, as in the example in Fig. ^] Moving a contour with filling ai past the pole at 
the origin changes its filling to a[ = ai + 1; this follows from the expression of the filling 
fraction as a contour integral (2.23), the residue theorem and the fact that the residue 
at X = switches its sign as the contour moves past it. As a result, the absolute density 
on the standard contour typically exceeds unity at its centre, as in Fig. 21, This and 

^Here and in the following, the cut contours that are determined through the requirement that 
p{x) dx e M, and which do not involve condensate cuts will be referred to as "standard contours" . 
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Figure 15: Magnification of tlie dotted area in Fig. 14 
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(a) Cut contours 




-0.04 -0.02 0.02 0.04 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 16: Cut contours and absolute densities of an example solution in region 
Ii2 (see Fig. |ll}. The fillings are ai = 0.016 (blue/thick, m = 1) and a2 = 0.03 
(red/thin, n2 = 2). 
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(a) Cut contours 
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(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 17: Another configuration in region I12. The cut contours and densities 
of Fig. 16 are shown with dashed lines. Cut one has a larger filling than in Fig. 16 
and attracts cut two (ai = 0.06, 02 = 0.03). 



(a) Cut contours 




-0.04 -0.02 0.02 0.04 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 18: A third configuration in region I12. The cut contours and densities 
of Fig. 16 are shown with dashed lines. Cut two has a larger filling than in Fig. 16] 
and attracts cut one (ai = 0.016, 02 = 0.06). 
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(a) Cut contours 



-0.10 -0.05 0.05 0.10 

(b) Absolute densities along the cuts vs. distance 
from the cuts' centres (measured along the con- 
tour) 



Figure 19: Standard cut contours and absolute densities of an example configu- 
ration in region II12 (see Fig. 11). The fillings fractions are ai = 0.11 (blue/thick, 
ni = 1) and 02 = 0.05 (red/thin, 712 = 2). This particular configuration lies be- 
low line G12 in the (a, E) plane. The physical version of the configuration has a 



condensate and is shown in Fig. 20 
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(a) Cut contours 



-0.10 -0.05 0.05 0.10 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the condensate and the cuts' contour) 



Figure 20: Physical version of the configuration shown in Fig. 19, The two 
common points of the standard contours are connected by a condensate. 



32 





(a) Cut contours 



-0.15 -0.10 -0.05 0.05 0.10 0.15 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 21: An example configuration in region II12 beyond line G12 (see Fig. 11 1. 
The filling fractions are ai = 0.032 (blue/thick, ni = 1) and 02 = 0.07 (red/thin, 
712 = 2) . The physical version of this solution has a condensate cut and is shown 



in Fig. 22 



the fact that the pole at the origin has the wrong sign in this configuration are strong 
indications that the physical realisation of such a solution is the one with a condensate 
cut, as in Fig. 22, 

Physically, there is no qualitative difference between configurations on either side of 
line G„^„+i, as their physical versions always contain a condensate cut and thus do not 
see the central part of the standard cut contours. Note that below line G„^„+i, as in 
the example in Fig. 19, the absolute density does not exceed unity on the standard cut 
contours. Hence the stability criterion (2.29) is obeyed by the standard contours and no 
condensate is necessary to fulfil (2.27). Yet a condensate version of the configuration is 
possible. Beyond line Gn,n+i a condensate is necessary for stability, and for continuity 
reasons one would assume that also below line G„ ,^+1 the physical version is the one with 
a condensate cut. This is another sign in favour of the assumption that if a condensate 
cut is possible, it is also realised physically. 

The other dashed lines H„^„+i and J„,n+i that further subdivide region Il^^n+i are 
discussed in Sec. ,,5.3, below. 

Region IIIn,n+i'- Two Disjoint Cuts, One with a Condensate. In this region, 
the contours of the two cuts are disjoint as in region I„,n+i, but the absolute density 
exceeds unity at the centre of cut two, i.e. the stability condition (2.29) is violated. An 
example configuration is shown in Fig. 23, Configurations in this region are still physical, 
because the same discussion as in Sec. 3.5 is applicable: A third, closed cut can be added 
in order to form a condensate at the centre of cut two. Since the absolute density at 
the centre of cut two exceeds unity, the excitation point with mode number n + 2 has 
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-0.06 -0.04 -0.02 0.02 0.04 0.06 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the condensate and the cuts' contours) 



(a) Cut contours 



Figure 22: Physical version of the solution shown in Fig. |21[ The part of cut 
one that encircles the origin is hidden in the condensate cut. 





(a) Cut contours 



-0.04 -0.02 0.02 0.04 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 23: An example solution in region III12 (see Fig. 11). The fillings are 
ai = 0.0038 (blue/thick, rii = 1) and 02 = 0.07 (red/thin, 77-2 = 2). The absolute 
density exceeds unity at the centre of cut two, hence the stability criterion (2.2S?j) 
is violated. The physical version of the configuration involves a condensate cut 
and is shown in Fig. 24, 
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(a) Central part of the 
cut contours 



-0.04 -0.02 0.02 0.04 

(b) Absolute density along cut two, the closed 
cut and the condensate vs. distance from the 
condensate's centre (measured along the con- 
tour) 



Figure 24: Physical version of the configuration shown in Fig. 23, A third, 
closed cut forms a condensate with cut one, as in the one-cut case discussed in 
Sec. |3.5| The closed cut starts and ends at the excitation point 1', which passed 
through cut two as the density at the centre of cut two reached unity. 



already passed through the contour of cut two and changed its mode number to n. The 
closed cut originates and ends at this excitation point. The resulting physical version of 
such a configuration is shown in Fig. 24j 

5.3 Lines 

In the following paragraphs, the features of configurations on the various lines shown in 
Fig. 11,12,13„14,15, are presented and appropriate examples are given. 

Line A.„: Only One Cut. This line represents configurations where only the cut 
with mode number n is present and has a filling a < Ocond (cf. Tab. 1), such that the 
absolute density at the cut's centre always stays below unity. Among all configurations 
with mode numbers n and n + 1 and given total filling, the state on this line is the one 
with the lowest energy. Adding a second cut with higher mode number while keeping 
the total filling constant obviously increases the energy. Conversely, adding a second 
cut with lower mode number while keeping the total filling constant decreases the total 
energy. 

Line B.^^: Only One Cut with a Condensate. This line is the continuation of line 
A„. Solutions on this line have only one cut with mode number n which has a filling 
ttcond < o. < ctcrit and wcrc discussed in Sec. 3, The absolute density at the centre of the 
cut's standard contour exceeds unity; hence a second, closed cut that starts and ends at 
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(a) 0.12 (b) 0.13 (c) 0.14 (d) 0.1487 (e) 0.1487 (f) 0.155 (g) 0.165 



Figure 25: A Series of configurations tliat passes through line Ci in Fig. 11 At 
all points of the series, the filling fraction of the small cut equals a^ = 0.04. The 
filling fraction ai of the large cut is indicated below each plot. Configurations 
(d) and (e) are the same in terms of mode numbers, fillings and branch points, 
but, as explained in the main text, this solution allows for two different contours 
for cut one. All solutions have a condensate in their physical version as displayed 
in Fig. i28| 



the excitation in — 1)' and forms a condensate with the central part of the cut is required 
for stability, as in the second example in Fig. 8j 

Line C^: Condensate with Two Tails. Line C„ separates region II„,„+i from region 



IL 



'l,n; 



which both consist of solutions with a condensate cut. Solutions on line Cr 



consist of a larger cut (cut one) with mode number n and a smaller cut (cut two), 
whose end points lie exactly on the standard contour of the larger cut. Therefore, the 
physical version of such a configuration consists of a condensate with one tail at each 
of its ends; the smaller cut is completely hidden in the condensate. Since the density 
of cut two at the cut's end points vanishes, the density of cut one at these points must 



equal — z, as explained below (2.14). Consequently, also the direction of cut one at these 
points is purely imaginary, i.e. vertical. Hence, the cut contour of the condensate version 
of this configuration and the density along it are smooth. Solutions on this line were 
investigated by Sutherland in |7|, refer to Figure 2 in that reference for a picture of the 
absolute density on the condensate and the tails. 

A series of configurations in which the branch points of cut two pass the contour of 
cut one is shown in Fig. 25. The sequence passes line Ci from left to right in Fig. 11, 
During the series, the filling fraction of cut two stays constant, while the filling fraction 
of cut one increases. At the beginning, the configuration is in region II12: The branch 
points of cut two lie to the left of cut one, and the mode number of cut two is ^2 = 2, 
while the mode number of cut one is rii = 1 (Fig. 25a, 25b, 25c). While the filling fraction 
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of cut one increases, cut two gets attracted by it and moves further right. During this 
process, the cuts increasingly bend towards each other and both become longer. At a 
certain value of the large cut's filling, the end points of the small cut lie exactly on the 
large cut's contour, the configuration has arrived at line Ci (Fig. 25d). One can see that 
the contour of cut one is indeed vertical at the end points of cut two. The filling of cut 
two stays constant during the series, but its length increases, hence the absolute density 
along its contour decreases. 

When the filling of cut one increases further, the small cut's branch points move on to 
its other side and the configuration reaches region II„_i,j (Fig. 25f,25g). While the end 
points of cut two pass the large cut, multiple quantities change: According to (4.13), the 
mode number of a cut is rii = iip^x^:), where x* is a branch point of the cut. Across a cut 



with mode number rii, p switches sign and shifts by 27mi (2.11). Therefore, a cut with 
mode number n2 whose branch points pass another cut with mode number ni changes its 
mode number to n2 = —n2 + 2ni. Specifically, if n2 = Ui + 1, as in the series in Fig. 25| 
then n'2 = ni — 1. Further, while cut two passes through the contour of cut one, the sign 
of the density (2.12) along cut two changes, as shown above in the discussion of region 
II„^„+i. This does not affect the density of roots along the cut, but it switches the sign 
of the filling: 02 — > 0^2 = —C(2- Thirdly, the contour of cut one and its filling change by a 
finite amount when cut two passes. The filling changes from ai to a'l = ai + 2^2, which 
is consistent with the fact that the total filling ai + 02, as a coefficient of the expansion of 



p at infinity (2.17), must vary continuously during the series. Also the total momentum 
m = riiai + 7120^2 varies continuously during the transition from region IIn,n+i to region 

J-J-n— l,n- 

For a qualitative picture of the densities on the cuts in the series, refer to Fig. '19b"26bl 
As a consequence of the change in the contour and the filling of cut one, in all solutions 
beyond line C„ (right of C„, in the (a, E) plane) the absolute density exceeds unity at the 
centre of cut one's standard contour. The standard contour of cut two is disjoint from 



cut one, as in Fig. |26[ In this form, the configuration does not allow for a condensate 
cut and is therefore unstable. However, there still is a stable version of such solutions: 
As discussed below (2.13), there are three possible starting directions for a cut at each 
branch point. Choosing a different direction for cut two, this cut's contour can be made 
to cross the contour of the first cut, and the unphysical density at the first cut's centre 
becomes hidden in a condensate. For the solution in Fig. 26, his type of configuration 
is shown in Fig. 27, These condensate-cut versions of the solutions beyond line C„ not 
only cure the instability at the centre of cut one, they also allow for a smooth transition 
from configurations on line C„ to solutions beyond line C„. This can be seen in Fig. 28^ 
where the physical version of the series of Fig. 25 is displayed. 

The changes that occur when the branch points of cut two pass the contour of cut one 
imply that configurations on line C^ have an ambiguity in the mode number of cut two 
and in the standard contour of cut one. This can be seen in the example in Fig. 25d,25e, 
In the limit where cut two approaches cut one from the left, the mode number of cut two 
is n2 = ni + 1, and the contour of cut one is bent to the left. In the other limit, where cut 
two approaches cut one from the right, the mode number of cut two is n2 = Tii — 1, and 
cut one is bent to the right. The fillings of the cuts behave accordingly, as described in 
the discussion of the series above. On line C„, the contour of cut one is unique between 
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(a) Cut contours 



-0.15 -0.10 -0.05 0.05 0.10 0.15 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 26: Standard cut contours and absolute densities of an example solu- 
tion in region II„_i^„, beyond line Ci (see Fig. 11). The fillings are ai = 0.25 
(blue/thick, ni = 1) and 02 = 0.04 (red/thin, n2 = 0). The absolute density 
exceeds unity at the centre of cut one, hence this configuration is unstable. The 
physical, stable form of this solution is shown in Fig. 27l 




(a) Cut contours 



-0.15 -0.10 -0.05 0.05 0.10 0.15 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the condensate and the cuts' contours) 



Figure 27: Physical version of the solution shown in Fig. 26. Cut two (red/thin) 
has a different starting direction than in Fig. 26 and hence a different contour, 
which allows for the formation of a condensate cut that hides the central part of 
cut one (blue/thick) and renders the configuration stable. 
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(a) 0.12 (b) 0.13 (c) 0.14 (d) 0.1487 (e) 0.1487 (f) 0.155 (g) 0.165 

Figure 28: Physical version of the series of solutions in Fig. 25. The configu- 
rations (f ) and (g) show the physical, smooth continuation of the series beyond 
line Ci. 



its end points and the branch points of cut two. The fact that the contour is ambiguous 
between the two branch points of cut two is consistent with the fact that the expansion 
of p has a square root term at the branch points x* (cf. (2.13)): 



p(x, +e)= p{x^) + A^e + 0{e) 



because this results in the expansion of the density (as follows from (J2.12)) 

p\{x^ + e) = ^ (p(x* + e) - vrni) = i (ni - n2 - -kA^/e) + 0{e) , 



(5.i: 



(5.2) 



which implies that the solution of the differential equation pi(x)dx G M has an ambigu- 
ous second derivative at the branch points x± and hence there are two choices for the 
continuation of the contour beyond the branch points x^,. However, this ambiguity has no 
physical significance, as the physical configuration with the condensate is non- ambiguous. 

Line D„^„_|_i: Two Tangential Cuts. This line separates region I„^„+i from region 
II„^n+i- The two cut contours of configurations on this line bend towards each other and 
touch each other in one point on the real axis. At the common point of the contours, 
the sum of the two densities equals —i. An example is shown in Fig. 29, Beyond this 
line, the formation of condensate cuts begins. 

Line E^^„_j_i: Cut Two Violating Stability. This line can be reached from region 
In,n+i by increasing the filling of cut two while keeping the filling of cut one small, such 
that the two cuts stay disjoint. Line E„^„+i is reached when the absolute density at the 
centre of cut one reaches unity. At this point, the excitation with mode number n + 2 
collides with the contour of cut two from the left. Further increasing the filling of cut 
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(a) Cut contours 



-0.05 0.05 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres (measured along 
the contour) 



Figure 29: An example configuration on line D12. The fillings are ai = 0.05833 

(blue/thick, m = 1) and 02 = 0.06 (red/thin, n2 = 2). The cuts touch at one 
point on the real axis, where the sum of their densities equals —i. 

two makes the excitation pass to the right of the contour and change its mode number 
to n. The absolute density on the standard contour of cut two now exceeds unity and 
a third, closed cut is required for stability; the configuration has reached region III„ „+i 
(cf. Fig. 24). 

Line Fn^n+i- Two Disjoint Cuts and a Cusp on Cut One. This line separates 
region III„„_|_i, where the standard cut contours of the two cuts are disjoint, from solu- 
tions in region Il^^n+i, in which cut one encircles the origin and consequently crosses the 
contour of cut two. A sequence of configurations that passes line F12 from region III12 



to region II12 is shown in Fig. 30, The transition happens when the excitation n', which 
carries a closed loop that forms a condensate with cut two, reaches the contour of cut 
one. At this point, the contour of the closed loop merges with cut one, whose standard 
contour encircles the origin afterwards. The transition point is marked by line F„^„+i. 
On this line, the density (2.9) on cut one decreases to zero at the merging point xq, as 
p{xo) = Tin. Around this point, the density expands to pi(xo + e)de ~ e^, which implies 
that the contour has a cusp with opening angle 90°. 



Line Gn,Ti-i-i: Condensate with Tails and a Virtual Cusp. This line in region 
II„ „+i separates solutions in which the standard contour of cut one closes on the right 
side of the origin from solutions in which it encircles the origin (cf. Fig. 19pT ). The 
qualitative change that happens to cut one's contour when crossing this line is the same 
as on line F„^„+i: The excitation point n' collides with the standard contour, making 
the density decay to zero at its centre. Unlike on line F„_„+i, here the process is less 
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(a) 0.0703 



(b) 0.07042 



(c) 0.0705 



Figure 30: A sequence of configurations that passes line F12. Cut one 
(blue/thick) has mode number ni = 1 and constant filling ai = 0.005, cut 
two (red/thin) has mode number n2 = 2 and its filling 02 is indicated below each 
subfigure. The excitation with mode number 1 is shown and carries a closed cut 
(green/thick) which forms a condensate cut at the centre of cut two. Configura- 
tion (a) lies in region III12, configuration (b) lies on line F12 and configuration (c) 
lies in region II12. The dashed contours all close around the origin. The energy 
increases during the sequence. 

significant because the central part of cut one is hidden in the condensate. An example 
configuration that lies on line G12 is displayed in Fig. _3]_, 

Line H„: One Cut with Vanishing Virtual Filling. On this line, the filling ai of 
cut one (with mode number n — 1) vanishes. However, as shown above in the discussion 
of region Il„,i_|_i, the standard contour of cut one encircles the origin and the filling on 
this contour equals one. Since most of cut one is hidden in the condensate cut, the 
vanishing filling of the physical cut contour does not have any further implications for 
the physical configuration. An example is given in Fig. [32| Because the filling of cut one 
vanishes, there is a one-cut solution with the same mode numbers and fillings for each 
solution on this line. These one-cut solutions are of the type shown in Fig. 8c and have a 
higher total energy than the corresponding solutions on line H„. Refer to the discussion 
of point Yn below for more details on the relation between these two types of solutions. 

Line Jn'- Excitation n Collides with the Condensate. These lines mark all con- 
figurations where the excitation point with mode number n is located on the condensate 
cut. When line J„ gets traversed, the excitation point with mode number n collides with 
the condensate from the left and passes to the right (for positive mode numbers), thereby 
changing its mode number to n — 2. But unlike an excitation point passing a standard 
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(a) Cut contours 



-0.05 0.05 

(b) Absolute densities along the cuts vs. dis- 
tance from the cuts' centres 



Figure 31: An example configuration on line Gi2. The fillings are ai = 0.0635 
(blue/thick, ni = 1) and 02 = 0.064 (red/thin, n2 = 2). The excitation with 
mode number 1 is shown and lies on the standard contour of cut one, which 
implies that the density (12.12) decays to zero at that point. 
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\ 
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(a) Physical version of the cut con- 
tours. Dashed lines show the hidden 
part of the standard contours. 




0.8 
0.6 
0.4 
0.2 



-0.04 -0.02 0.02 0.04 

(b) Absolute densities along the condensate and 
the tails vs. distance from the condensate's cen- 
tre (measured along the contour). 



Figure 32: An example configuration on line H2. The fillings are ai = 
(blue/thick, ni = 1) and 02 = 0.085 (red/thin, 77-2 = 2). The parts of the 
standard contours that are hidden in the condensate cut are shown as dashed 
lines. Since the standard contour of cut one encircles the origin, the integral of 
the absolute density along it is not zero, but unity. 
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cut, this process does not create an additional, closed cut (as on line E„^„+i); in fact, the 
passage of the excitation point has no effect on the configuration at all. 

5.4 Special Points 

As a result of the discussion of the various regions and lines in the previous sections, three 
types of special points emerge, labelled by X„, Y^ and Z„_„+i in the Fig. Il|12[13[14p5 
These points will be discussed in the following. 

Point Xrii Excitation Point Meets Single Cut. The points X„ mark the one-cut 
solutions with mode number n whose absolute density equals unity at the cut's centre, 
which means that the excitation point with mode number n + 1 is located on the cut's 
contour. The filling fractions of these solutions for n < 10 can be found in Tab. 1, The 
point X„ connects line A„, on which the absolute density on the single cut's centre falls 
below unity, with line B„, on which the absolute density at the centre of the standard 
contour exceeds unity and a condensate begins to form with the help of a second, closed 
cut. Also the lines C„, Y)n,n+i and E„_i „ begin at the point X„ (the last one only if 
n > 1). Among all configurations on line A„, the one at this point has the highest energy 
(cf. Tab.|lJ). 

Point Y„: The Two Excitation Points with the Same Mode Number Collide. 

Point Y„ marks the end of line B„ and thus represents a one-cut solution with a conden- 
sate. On line B„, the two excitation points {n — l)' and (^ — 1) both lie on the real line, as 
in Fig. |8bj At the point Y„, these two excitation points collide. Analytically continuing 
line B„ beyond point Y„ lets the branch points diverge into the complex plane and leads 
to one-cut solutions of the type shown in Fig. |8c. This analytic continuation of line B„ 



shall be called line B^. As discussed in Sec. 3.5, solutions on this line are limiting cases 



of three-cut solutions. Hence line B^ is not part of the moduli space of two-cut solutions 
(or of its boundary). 

The only way to continue beyond point Y„ while staying in the regime of consecutive 
mode numbers is to enter either region II„_i„ or region III„_i,^. In both cases, the 
excitation points {n — 1)' and (n — 1) separate again. The former stays on the real line, 
while the latter splits into two square root branch points that diverge into the complex 
plane. Upon entering region III„_i^„, the excitation point {n — 1)' continues to carry a 
closed cut as on line B„ while the new branch points form the ends of a second standard 
cut, as in Fig. 30a, When the configuration enters region II„_i .^, the afore closed cut 
opens up at its cusp and becomes a standard cut that ends on the new pair of branch 
points, as in Fig. 30c, Either way, a second standard cut emerges, which can have a 
positive, negative or vanishing filling. Configurations with vanishing filling belong to 
line H„ and lie in region II„_i „. Since the closed cut in the solutions on line B„ also has 
zero filling, line H„ seems to be the most natural continuation of line B„. 

As noted at the end of Sec. 3, solutions on line B'^ can be transformed into the 
corresponding solutions on line Hn by separating the two segments of the closed cut 



in Fig. 8c, In doing so, the two excitation points in the complex plane split into two 



43 



n ^ 

47r^ 



-0.001 



-0.002 



B2 : 


Y2 1 B^ i 


i 


hN^ 


\ 


X 


1 1 


1 1 1 1 1 1 m > 



-0.005 







0.005 



0.01 



0.015 a — ttYa 



Figure 33: Deviation of the energy on line H2 from the energy on hne Bj 
(green/dashed), which is the analytic continuation of line B2 (purple/solid). 
Shown are the energies of numerically obtained solutions on line H2 (black dots), 
the corresponding function £'h„ — -Eb' (purple/solid) and the quadratic term of 



(5.3) (purple/dashed). The energy deviation is measured between solutions with 
the same total filling a and shows that a second order phase transition occurs 
upon the passage from line B„ to line H„. 



square-root branch points each, which results in a three cut solution. The filling of 
the third, vertical cut can then be decreased to zero while the filling of the second 
cut increases to zero, such that the total filling fraction a stays constant. During this 
process, the energy decreases. Comparing the energy of solutions on line H„ to the 
energy of the corresponding one-cut solutions on line B^ (Fig. 33) shows that a second 
order phase transition occurs as a configuration follows line B„ and continues on line H„. 
For example, the deviation of the energy on line H„ from the energy on line B^ has the 
following expansion in the total filling at point Y„: 



(^H„-^B')(«) 



-K n {a- a;Y,J + 0{{a - ay^ 



(5.3) 



The prefactor of the quadratic deviation, including its dependency on n, has been guessed 
on the basis of numerical data which, for the first few values of ra, matches this form with 
high accuracy. On a technical level, this phase transition is related to the one discussed 
by Douglas and Kazakov in |29|: It involves a transition between the same classes of 
algebraic curves and it is also related to a density threshold, see App. '^ for a technical 
review. The authors study the large- iV limit of the partition function of continuous pure 
Yang-Mills theory on the two-sphere and find a phase transition with respect to the area 
of the sphere.^ 



''Wc thank V. Kazakov for remarking this point. 
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Point Zn,n+i'- Cusp on Cut One Meets Cut Two. At this point, the hnes F„ ,^+i 
and Gn,n+i meet: Cut one has a cusp at the centre of its contour whose tip hes on the 
contour of cut two. Since the tip of the cusp marks the excitation point with mode 
number n + 2, this imphes that the absolute density of cut two at this point equals unity 
and that lines En^n+i and J n+2 also end on point Zn^n+i- The two cuts touch each other 
only in one point, thus also line D„_„+i ends here. 

5.5 Global Structure of the Moduli Space 

As a consequence of the discussion about the hnes C„ in Sec. 5.3, all moduli spaces of two 
consecutive mode numbers are analytically connected through the lines C^. Pictures of 
the resulting global space are shown in Fig. [34p5| In the latter, also the lines of unstable 
one-cut solutions B^ are shown for comparison, while in the former these lines are the 
same as the lines H„. In particular, the global connectedness of the space shows that 
mode numbers lose their intuitive meaning at large filling fractions: A small change in 
the number of excitations (filling fractions) can alter the mode numbers, hence they can 
no longer be interpreted as the number of periods of a coherent state on the spin chain. 

Note that the case rii = (cut one), n2 = 1 (cut two) qualitatively differs from 
the other cases with consecutive mode numbers. The excitation point ra^, = always 
lies at a; = cxD and hence cannot be excited, only the excitation point n^, = 0' can. 
Consequently, line Bi cannot be crossed, and cut one always has a negative filling and 
its standard contour winds around the origin. That means that solutions with mode 
numbers ni = 0, n2 = 1 always have a condensate and regions Iqi and IIIqi do not exist. 

The one-cut lines A„ and B„, n > 1, form branch cuts of the global moduli space; 
across them, the derivatives of physical quantities such as the energy and the momentum 
are discontinuous. The points Y„, n > 1 are the corresponding branch points. 



As mentioned at the end of Sec. 3^, at point Yi the excitation point n* = 0' reaches 
X = oo. At that point, the filling fraction of cut two equals 0:2 = 1/2 and cut one extends 
along the whole imaginary line while the branch points of cut two are x± = ±i/27r. The 
absolute density of this solution is shown in Fig. 36. When the filling of cut one gets 
reduced to finite negative values while 0^2 stays at 1/2, the condensate and the tails of 
the configuration remain on the imaginary line, but their length decreases. Solutions of 
this type constitute line Co; an example is shown in Fig. 37, All solutions on this line 
have a momentum of P = 27r(niai -|- n2a2) = vr. 

As mentioned, a consequence of the fact that the phase space of consecutive mode 
numbers is connected is that the mode numbers can change by a continuous variation of 



the fillings. In particular, if one follows a straight horizontal line in Fig. |34jfrom the left 
to the right, the branch points cycle around each other while the mode numbers decrease 
until the configuration reaches line Co and both pairs of branch points and the contour 
lie on the imaginary line. A schematic sketch of this process is shown in Fig. 38. The 
sequence can even be continued further by following it in reverse, but in the other (left) 
half of the complex plane: The mode numbers become negative, but all other physical 
quantities are the same as in the mirror solution with positive mode numbers. This 
means that there exist two copies of the space in Fig. 34, one for positive and one for 
negative mode numbers, and the two parts are connected through line Cq. 
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Figure 34: Global structure of the moduli space for consecutive mode numbers, 
shown in the (P, P/a) plane. 
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Figure 35: Global structure of the moduli space for consecutive mode numbers, 
shown in the (P, E) plane. 
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Figure 36: Absolute density along the condensate and the two tails of the 
solution at point Yi. Cut one (ni = 0, ai = 0, red/thin) extends along the 
whole imaginary line, while the branch points of cut two (n2 = 1, 02 = 1/2) are 
x± = ±i/2-K. The absolute density on the tails falls off as l/|x|. 
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Figure 37: An example solution on line Co- The condensate and the tails lie 
on the imaginary line. The standard contours are shown as dashed lines. The 
filling fractions are ai = —0.14 (red/thin, ni = 0, x± = ±0.2726^) and 0:2 = 1/2 
(blue/thick, n2 = 1, x± = 0.146434z). Cut one is completely hidden in the 
condensate. 
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Figure 38: A sequence of configurations that evolves through the global moduli 
space of consecutive mode numbers. During the sequence, the branch cuts cycle 
around each other. The mode numbers are indicated next to each branch point. 
Each time a branch point passes a contour, its mode number decreases by two. 
The total filling fraction a increases during the sequence; the individual filling 
fractions are indicated below each subfigure. The sequence roughly follows a 
straight horizontal line in Fig. 35 at E/Att'^ « 0.3. The final configuration lies 
on line Cq, i.e. all branch points lie on the imaginary line. The size of the cuts 
increases from left to right, the rightmost configuration is about 15 times larger 
than the leftmost. 



Also the lines B^, n > of unstable one-cut solutions (as shown in Fig. 8c) can 
be followed until the branch points reach the imaginary axis, similar to the sequence 
of two-cut solutions that approaches line Cq- At this point, also the fluctuation points 
n^: = {n — 1) and n* = (n — 1)' lie on the imaginary line and the fllling of the branch 
cut reaches a = 1/2, as follows from (3.5). But, according to (3.3), this implies that 
P = Tin > TT. Since the maximal physical momentum on a lattice is P = vr, the physical 
momentum is P — 27r for P > vr. The maximal physical momentum P = vr is reached by 
the configurations on line B'^ long before the branch points arrive at the imaginary axis. 
Reaching the regime of such high momenta is prevented by the phase transition discussed 
at the end of Sec. "5.4", which makes the one-cut solution transform into a two-cut solution 
at the critical value of the filling Ocrit = «¥„• Physically, nothing peculiar would happen 
when the momentum would increase to values larger than vr, but it is interesting to see 
that the two-cut solutions by themselves implement the restriction to momenta \P\ < it, 
as can be seen in the pictures of the phase space in Fig. 34p5 
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5.6 Stability and Gauge Theory States 

The expositions of this section show that the moduh space of two-cut solutions has a 
very rich structure, even though only solutions with consecutive mode numbers have been 
considered. Most remarkably, there seem to be no unstable two-cut solutions: For all 
spectral curves with consecutive mode numbers there is a configuration of standard and 
condensate cuts that generate the curve and that do not violate the stability condition. 
Here, by stability condition is meant that the absolute density never exceeds unity at 
any point on the cut contour whose tangent is vertical: 

|p(a;o)|<l forp(a;o)~^ (5.4) 

Such critical point are always hidden in a condensate cut. Thus it may be conjectured 
that all classically admissible two-cut solutions are solutions of the Bethe equations in 
the thermodynamic limit. In fact, numerical tests show strong evidence in favour of 
this conjecture: For most of the cases discussed in the previous section, numerically 
exact Bethe root distributions that approximate the cut contours have been obtained, 
see Sec. 

A remark about suitable gauge theory states is appropriate at this point. Only spin 
chain states with a momentum P G 27rZ correspond to gauge theory states. None of the 
solutions studied in this chapter satisfies this condition: All two-cut configurations with 
consecutive mode numbers have a non-zero momentum < P < vr. Nevertheless, all 
quahtative results should generalise to valid gauge theory states. Namely, configurations 
whose cut structure is point symmetric about the origin manifestly have a vanishing 
momentum P.^ It is reasonable to expect that such point symmetric solutions with 
mode numbers {n, n + 1) and (— n, —n — 1) and with symmetric fillings exhibit the same 
qualitative features as the solutions investigated here (such as condensate formation and 
connectedness of the moduli space). Some point-symmetric configurations of two-cut 
and degenerate multi-cut solutions are studied in the next section. 

6 Moduli Space of Non-Consecutive Mode Numbers 

All qualitative features of the general moduli space with arbitrary mode numbers can 
already be seen in the case of consecutive mode numbers, which was studied in the 
previous section. In general, cuts with non-consecutive mode numbers lie further away 
from each other and hence influence each other less. Some examples will be studied in 
this section. 

6.1 Mode Numbers ni = 1, n2 = 3 

First, the comparatively simple space of solutions with mode numbers ni = 1 (cut one). 



71.2 = 3 (cut two) shall be exposed. It is shown in the (a,P/2Tca) plane in Fig. 39 



In region I, the fillings of both cuts are small, their contours are disjoint and their 



^Unless there is a condensate cut on the imaginary line that passes through the origin. In that case, 
the niomcntuni is maximal, P = tt. 
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Figure 39: The space of solutions witli mode numbers ni = 1, n2 = 3, sliown 
in tlie (a,P/27ra) plane. Some example solutions in different regions are shown. 
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absolute densities are bound by unity, hence there are no condensates. Line Di marks 
the configurations in which the density on cut one equals —i. Beyond this line lies region 
III, whose configurations have a condensate on cut one and an associated closed loop-cut 
with mode number = 2'. Similarly, configurations on line D2 have an absolute density 
of —i at the centre of cut two; beyond this line lies region II3, where cut two has a 
condensate core and an associated closed loop-cut with mode number 2 = 4'. In region 
II13, both cuts have a condensate core and associated closed loop-cuts. Finally, line C24 
marks configurations in which the excitation points with mode numbers n = 2 and n = A 
have collided. Configurations beyond this line are unstable in the sense that they are 
degenerate cases of solutions with higher genus (higher number of cuts). 

Regions IIi and II13 can be further subdivided, which is signified by the dashed 
and dotted lines in Fig. [39J To the left of these lines (at low total filling a), the two 
cuts with their condensates and loop-cuts are disjoint. On the green, dashed fine, the 
excitation point with mode number 2' = 4 (which is the cusp of the loop-cut of cut two) 
collides with the condensate core of cut one. Upon further increasing the total filling, the 
condensate core of cut two (or its standard contour) collides with the condensate core of 
cut one; this event is marked by the orange, dotted line. At even higher total filling, the 
branch points of cut two pass through the condensate core of cut one (orange, dashed 
line), such that the whole contour of cut two (with its associated loop-cut) is encircled 
by the loop-cut of cut one. Nothing special happens to the condensate core of cut one 
during the passage of cut two - this is sensible, as the condensate core is invisible in 
the spectral curve dp, which contains all the information about the solution. When cut 
two has completely passed through the condensate core of cut one, its mode number has 
changed to n'g = n2 — 2 = 1. 

When the total filling gets increased up to a = 1/2 (line E_i), the excitation point 
with mode number 2 = 0', which marks the cusp of the closed loop-cut that is connected 
to cut one, reaches x = cxd. This results in a configuration similar to the one at the 
point Yi (see Fig. 36): The contour of the closed cut passes through the branch points 
of cut one, such that cut one is completely hidden in a condensate that has two vertical, 
infinitely long tails. Consequently, the mode number of cut one is ambiguous: It is either 
rii = 1 or rii = —1. In addition, there is cut two, whose mode number is n2 = 1. It has 
a standard contour (below line D3), or a standard contour with a condensate core and a 
closed loop-cut (above line D3). At the upper end of line E_i lies point V, at which also 
the excitation point that marks the cusp of cut two's closed loop-cut has reached x = 00 
and the configuration is completely symmetric about the origin (Fig. 40a). Consequently, 
the momentum vanishes at this point, P = 27r, and the solution represents a valid gauge- 
theory state. Fig. |40b| shows another cut configuration that realises the solution at point 
V. It can be constructed by taking another possible standard contour for each cut and 
replacing their common centre with a double condensate, that is a condensate with root 
density \p\ = 2. Bethe root distributions that show this type of pattern have been 
found before |9,8|. Behind point V lies line Ei. Configurations on this line are "mirror 
solutions" of the ones on line E_i: Cut three has mode number 1/ — 1 and consists 
of a condensate with two vertical tails that extend to infinity, while cut one has mode 
number —1 and consists of a standard contour with or without a condensate core and 
an associated closed loop-cut. 
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(a) (b) 

Figure 40: Configuration at point V in Fig. 39, This solution is symmetric 
about the origin and hence represents a valid gauge-theory state (The momentum 
equals P = 27r). Figure (a) shows the configuration that is reached when coming 
from region II13 or from line E_i in Fig. ^9^, There is a dual configuration of cuts, 
shown in '(b) , that generates the same solution. It has a "double-condensate" with 
root density |/j| = 2 at it's centre. 



6.2 Mode Numbers rii = 1, 712 = — 1 

The moduli space of solutions with mode numbers rti = 1 (cut one), n2 = — 1 (cut two) 
is shown in Fig. [41p2j in the (a,P/2'7rQ;) and in the (a,-E) plane. Regions, lines and 
points with the same labels but subscripts 1 and —1 fall together in Fig. 42, since their 
configurations only differ by a reflection about the imaginary line and hence have the 
same energy. 

As in previous examples, the lines A±i and B-ti represent solutions with a single 
cut with mode number n, solutions on the lines A±i have a standard contour without a 
condensate core, while solutions on the lines B-ti have a condensate core and a closed 
loop-cut. In the regions I-ti and II±i, both cuts with mode numbers rti = 1 and n2 = —1 
are present. Line D±i separates solutions that only have a standard contour (region I-ti) 
from solutions that have a condensate core with a closed loop-cut (region II±i). The 
closed loop-cut has mode number n = for both cuts. Hence on line 'D±i, the absolute 
density at the centre of the cut with mode number ±1 equals \p\ = 1. 

In the regions Ii and IIi, the flUing fraction of cut one is larger than the fllling fraction 
of cut two, ai > a2- In the regions I_i and II_i, it is the other way round, ai < 02. These 
regions are separated by the lines S and Sc, which consist of symmetric solutions with 
«! = a2- Such solutions have vanishing momentum P and hence represent valid gauge 
theory states, they were studied in |27|. Obviously, solutions on line S have plain standard 
contours, while solutions on line Sc have condensate cores in both contours. Point Xg 
marks the symmetric solution with maximal fllling that does not have condensates. The 
lines E-ti and point V are the same as the ones discussed above in Sec. 6.1, hence the 
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Figure 41: The space of configurations with mode numbers ni = 1, n2 = —1, 
shown in the (a, P/27ra) plane. Configurations in regions I±i have plain standard 
contours, while configurations in regions Il-ti have condensates on either or both 
of the contours. 




0.5 « 



Figure 42: The space of configurations with mode numbers ni = 1, n2 = —1, 
shown in the (a, E) plane. Configurations that only differ by a reflection about 
the imaginary line have the same energy, hence the upper and the lower lines 
and regions of Fig. 41 lie on top of each other in this representation. 
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diagrams in Fig. 39141 are connected. The points Y-ti actually both represent the same 



solution, namely the one shown in Fig. 36 whose cut configuration lies entirely on the 
imaginary line. 

7 Numerical Solution of the Bethe Equation 

In this section we develop a method for obtaining numerically exact distributions of Bethe 
roots for the Heisenberg spin chain. The main goal of this section is to test our predictions 
from the analysis of the one-cut and two-cut spectral curves in Sec. 3,5,6, Thus we are 
especially interested in distributions which can be compared with the spectral curve, i.e. 
solutions with large numbers of roots that align on a set of contours. 

We examine in detail some subtle points of the previous analysis, such as stabil- 
ity (Sec. 2.7), condensate formation (Sec. 3.5) and phase transitions between spectral 



curves with different topologies (Sec. |5.4 ). We will see that condensate cores with roots 
separated by Au ^ i indeed appear in all cases considered. We will also confirm that 
the phase transition from one-cut to two-cut solutions indeed happens. Furthermore, in 
Sec. 7.3 we will analyse small deviations from the thermodynamic limit. In particular 
we will numerically analyse the leading finite-size corrections and show that they do not 
exhibit any non-analyticity due to the formation of the condensate. 

The numerical configurations shown in this section are approximate solutions to the 
discrete Bethe equations (1.1) with finite length L. It is astonishing to observe how well 
the predictions from the thermodynamic limit describe the actual distributions of Bethe 
roots. 

7.1 Description of the Numerical Method 

Here we describe the numerical method used to solve the Bethe equations (1.1) in detail. 
In order to find roots of a set of equations numerically, it is important to specify a good 
set of starting points that are sufficiently close to the exact solution. We will focus on 
this problem first and then describe a simple programme for finding numerical solutions. 

Starting Points. In Sec. '2.3' we introduced the mode number associated to a branch 
cut in the thermodynamic limit. We will need to generalise the notion of mode numbers 
to the quantum case. The mode numbers naturally appear in the quantum theory when 



one takes the logarithm of (ll.l): 



tLlog'^ + 2nn, = Vzlog ^^ ~ ^^^ + ^. , k = l,...,M. (7.1) 

The integers Uk represent the 27rz-ambiguity of the logarithm, in the thermodynamic 



limit they exactly become the mode numbers introduced in equation (2.11). This can be 



seen by taking the large-L limit of (7.1), as in [7|[8|,[18|. The solution to (7.1) is uniquely 



specified by a set of integers n^ (at least for large L and small M). 
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Figure 43: Approximation of a solution to the Bethe equations (1.1) by zeros 
of the Hermite polynomial for M = 16 roots with mode number n = \ and 
chain lengths L = 400, L = 1600. Squares represent points given by Uk = 
L/(27rn) + izfc\/2^/(2vm), where Z]^ are zeros of Hermite polynomial of degree M, 
while circular dots represent numerically exact Bethe roots. This approximation 
provides us with a good set of starting points. 



In order to find the solution numerically by Newton's method, one has to specify 
starting points for the search. For generic parameters it is very difficult to give suitable 
starting points that lead to convergent iterations. However, as we will see it is pretty 
simple to find good starting points in the limit of small filling fractions L ^ M where 
the Bethe equations (jYTj) can be solved approximately |9|: 

To the leading order we can drop the r.h.s. of (7.1) and, assuming that m^ ~ L (which 
is true in the case we are interested in), we get that Uk — L/27mfc. We conclude that 
Bethe roots with the same mode number Uk are close to each other and separated by 
Am ~ L from other groups of roots with different mode numbers. For resolving the 
positions of the Bethe roots with the same mode number n, we can neglect terms in the 
r.h.s. of (7.1) containing roots with different mode numbers and use the ansatz 



Uk 



Inn \ 



Assuming that n ~ 1 and z^ ~ 1, we expand (7.1) and get |9 



^k Zj 



(7.2) 



(7.3) 



which is the equation for the zeros of the Hermite polynomial, HM(-2fc) = 0. Numerically 
exact values of Bethe roots are plotted against their large-L approximation (|7.2[) in 
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Fig. 43, and one clearly sees that the relative mismatch is indeed smaller for larger 
L. It is possible to continue the l/yL expansion in a systematic way to get a better 
description of Bethe roots for large L and mode number n fixed. For our goal the leading 
approximation of the Bethe roots by zeros of Hermite polynomials provides us with a 
good set of starting points for sufficiently large L. 

Iteration Procedure. Having starting points at hand it is almost immediate to con- 
struct a numerical solution of the Bethe equations using Newton's method. However, we 
know the starting points only for large L. To overcome this difficulty we can view the 
Bethe root distributions as functions of the length L while keeping the magnon number 
M fixed. We can solve the Bethe equations numerically, gradually varying L from a 
sufficiently large value L :» M to the desired value with L = 0{M). In each step we use 
the numerical solution obtained in the previous step as new starting values for smaller L. 
In this way we obtain a series of configurations whose filling a = M/L slowly increases. 
Some important examples of such series will be considered in the next subsection. 

For various reasons it is favourable to use the Bethe equations in product form pTl| 
and not in logarithmic form (7.1) for the Newton iteration. Using this method, the 
mode numbers enter into the calculation only through the starting points. Since we 
vary L slowly, at each step we inherit the information about the mode numbers from 
the previous step, and hence the resulting configuration corresponds to the analytically 
continued curve from the regime of large L. 

Note that for large L and fixed M the roots are separated by Au ~ vZ and thus the 
stability condition is satisfied, \p\ ~ 1/An ^ 1. However, when we decrease L, at some 
point the separation between Bethe roots approaches Au = i. The Bethe equations ^.Ij) 
are singular at this point which leads to a numerical instability at small L.^ Therefore 
the step size for L must be chosen sufficiently small in this regime. For practical purposes 
one can even choose L to take real intermediate values. We expect the above numerical 
procedure to lead to a reasonably unambiguous final configuration if L is changed in 
sufficiently small steps. 

Sample Programme. Let us give a basic realisation of above method in Mathematica. 
For simplicity the programme presented in Tab. 2 is restricted to the one-cut case. 
To generalise it, it is enough to change the set of starting points uO according to the 
discussion in the previous section. We tested this programme on Mathematica versions 
5.2 and 6.0, but it should be compatible with earlier versions as well. 

This programme can be optimised to work much faster, especially for the case of a 
large number of Bethe roots: A major improvement can be achieved by replacing the 
built-in FindRoot function, which is too general and inefficient, by a simple implemen- 
tation of Newton's method. Furthermore the step size for L can be made adaptive. 
One can even decrease L in small steps during each Newton iteration and thus avoid an 
excessively precise calculation which is merely used for new starting values. 



^In the logarithmic Bethe equations (7.1) the branch cut of the logarithm furthermore introduces a 
discontinuity close to the starting values. The Newton method is thus likely not to converge and the 
product form (1.1) is clearly preferable in this case. 
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(♦define mode number, number of roots, initial length and finial length*) 
n = 1; M = 12; L = 100; Lfinal = 30; 

(♦starting points*) 
uO = (L/(2 Pi) + 1 Sqrt[2 L]/(2 Pi n) zk / . NSolve [HermiteH[M, zk] == 0, zk] ) 

(♦solve the Bethe equations while slowly decreasing the length L*) 
Do[uO = Table[u[k], {k, M}] /. FindRoot [Table [((u[k] + l/2)/(u[k] - 1/2))"L == 
-Product [(u[k] - u[j] + l)/(u[k] - u[j] - 1), {j, M}] , {k, M}] , 
Table[{u[k], uO[[k]]}, {k, M}] , WorkingPrecision -> 50, AccuracyGoal -> 20]; 
L— // Print; 
, {L - Lfinal + 1}] 

(♦plot the result*) 
ListPlot [Transpose [{uO//Re,uO//lm}] , PlotRange -> All, AspectRatio -> Automatic] 

Table 2: Mathematica programme for the generation of a finite-length numerical 
Bethe root distribution corresponding to a one-cut spectral curve. 

7.2 Some Series of Numerical Solutions 

In this section we will exploit numerical solutions of the Bethe equations obtained with 
our programme to test the predictions from the thermodynamic limit. We begin with a 
study of the simplest one-cut solution, subsequently proceed to more complex configu- 
rations with two cuts and consecutive mode numbers, and finish with a two-cut solution 
with non-consecutive mode numbers. In the last configuration, one of the two cuts 
transforms into a two-cut complex in a phase transition, the result being a three-cut 
solution. 

One-Cut Solution, n = 1. In the thermodynamic limit this solution is discussed in 
Sec. [3j An interesting phenomenon which we expect to arise for this configuration is the 
appearance of a closed loop-cut for filling fractions above ctcond (see Fig. 8b). 

The spectral curve of the one-cut solution can be discussed very explicitly, in partic- 
ular the quasi-momentum is given by a simple algebraic formula (3.1) and one can find 



the shape of the physical branch cut quite easily. In the first figure of Fig. j44jthe contour 
is shown together with the numerical solution for M = 60 roots and length L = 2925. 
We almost do not see any deviation of the roots from the asymptotical curve. On the 
next picture we considerably increased the filling fraction to M/L = 60/315 ~ 0.19, and 
one can see that the distance between roots close to the real axis is just slightly above i. 
This means that the filling a is very close to acond and the stability condition is almost 
violated. According to the discussion in Sec. j3.5| this implies that the fluctuation with 
mode number 2 is very close to the cut. In the figure, we indicate it by a red diamond. 

When we continue increasing the filling fraction to values above Ocond, the fluctuation 
point passes through the cut contour and produces a loop cut. At the intersection point 
of the loop cut with the original standard cut we expect the condensate to be formed. 
As is demonstrated in Fig. 44, this is indeed the case. 

We can almost reach the maximal filling a = \ numerically. As discussed in Sec. "5.5l, 
this would correspond to a configuration with a condensate centred at the origin together 
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Figure 44: One-cut solutions with mode number n = 1 and lengths L = 
2925,315,160,130. Shown are the cut contours of the spectral curve together 
with numerical solutions. When the maximum density along the cut exceeds 
\p\ = i, the fluctuation point with mode number n = 2 touches the cut. When 
we continue decreasing L, a condensate develops in the middle of the original cut 
and two extra tails appear. In fact in the scaling limit L ~ M the two tails to 



the right are touching at the point of the virtual excitation, cf. Sec. 3.5 



with two tails that extend along the imaginary axis towards the fluctuation point, which 
has diverged to infinity. 

Finally, we notice that the roots belonging to the two tails on the right of the conden- 
sate do not seem to align on the contour predicted by the spectral curve (red/dashed). 
In fact this is an artefact of the finiteness of the number of roots. One can observe that 
the situation becomes better when the number of roots M is increased (while the filling 
fraction M/L is held fixed). The finite-size effects become very relevant at the tip of the 
loop because the density drops to zero at the fluctuation point forming the tip. Close 
to a standard branch point the roots align on the contour more accurately, because the 
density decreases only as the square root of the distance from the branch point, whereas 
at the cusp of the loop the density is proportional to the distance itself and hence the 
deviation from the spectral curve is more severe. 



One-Cut Solution, n = 2. For larger mode numbers, the one-cut solution behaves 
similarly to the n = 1 case when the filling fraction is increased. However, we expect 
that a configuration with n > 1 exhibits a new type of behaviour when its filling fraction 
reaches aait given by ([3^). As discussed in Sec. 3.5, we expect that a phase transition 
from one cut to two cuts occurs. 

In Fig. "45", we compare the distributions of numerical Bethe roots with mode number 
n = 2 for different filling fractions. We see that the condensate forms when the nearest 
fluctuation from the left crosses the cut contour at a = Ocond, which is in full agreement 
with the prediction of Tab. 1. In contrast to the n = 1 case, there is also a fluctuation 
point to the right of the cut (with mode number n = 1), even at small fillings a. When 
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Figure 45: The one-cut configuration with mode number n = 2 transforms 
into a two-cut configuration after the collision of the two fluctuation points with 
mode numbers n = 1, 3. 
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Figure 46: Shown is the same diagram as in Fig. 33, In addition, the energies of 
a series of numerical solutions to the discrete Bethe equations that start as a one- 
cut solution with mode number n = 2 (red/solid) are displayed. The graph shows 
that the one-cut solution (line B2) indeed transforms into a two-cut solution (line 
H2) when the two fluctuation points collide (point Y2). The deviation from the 
exact thermodynamic lines is due to the finite length and number of roots in the 
numerical solutions. 
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Figure 47: Two-cut solutions witli consecutive mode numbers n = 1 (riglit, 
blue) and n = 2 (left, red) for different values of the length L. The intersection 
points of the two standard contours are connected by a vertical condensate cut. 



we increase a, the fluctuation point with mode number n = 3, which initially was on the 
left of the standard cut, passes through the cut and approaches the fluctuation point with 
mode number n = 1 (second flgure in Fig. 45). Beyond the point of collision there are 
two possibilities - the conflguration could continue as a one-cut solution (with a loop cut 
with two cusps that could naturally be decomposed into two standard cuts), or it could 
transform into a two-cut solution with lower energy as explained in Sec. |3.5[ We found 
that in our study always the second possibility is realised, our numerical method prefers 
the two-cut conflguration with lower energy. It seems that this conflguration is numeri- 
cally more stable. In particular we see that in the last two pictures of Fig. 45 the roots 
perfectly follow the two standard contours of the two-cut spectral curve. Accordingly, 
Fig. 46 shows that the energy of the numerical conflguration after the point of colli- 
sion matches the energy of the two-cut conflguration (up to deviations due to flnite-size 
effects), which clearly indicates that the phase transition is indeed taking place. 



Two-Cut Solution, n = 1,2. The spectral curve for two-cut solutions was con- 
structed in Sec. 4 and is less simple than the one-cut curve. The moduli space of two-cut 
conflgurations is also much more complicated, for the case of consecutive mode numbers 
it is described in detail in Sec. 5, Numerically, however, it is not difficult to construct 
a solution with an arbitrary number of cuts. We mainly want to show here that, in 
accordance with the prediction of Sec. 5, an intersection of the two cuts leads to the 
formation of a condensate that connects the two intersection points. 

In Fig. 47 we compare numerical two-cut solutions with mode numbers n = 1, 2 with 
the predictions from the spectral curve. In the first two pictures the two cuts are disjoint 
and hence we are in region I12 of the phase diagram shown in Fig. 11. The last three 
configurations belong to region II12, and we see that indeed a condensate with roots 
separated by approximately i connects the intersection points of the two contours. 



61 



L=35()M,=3()M,=1() 



L=22()M|=3()M2=1() L=140M|=30M,=10 L=114M|=30M,= 10 



L=89M|=30M2=10 



20 



.1=4! 




20 30 40 50 



15 


• 




• 


10 




5 









-5 




-10 






• [ 


-15 


• 



15 


• 




10 




".• 


5 













-5 






-10 


•' 




-15 







10 15 20 25 30 6 8 10 12 14 16 2 4 6 8 10 12 2 4 S 8 10 12 14 



Figure 48: Two-cut solutions with non-consecutive mode numbers n = 1 (blue) 
and n = 3 (red) for different values of the length L. Both cuts develop con- 
densates before they intersect each other. We see that the passage of one cut 
through the condensate core of the other does not affect the relative positions of 
the roots inside each cut. 



On the third picture, the right (blue) n = 1 cut does not encircle the origin of the 
complex plane, but it has a cusp at the point of intersection with the real axis. This 
means that this configuration lies exactly on the line G12. In the last configuration shown 
in Fig. |47l the branch points of the n = 2 cut (red) are almost touching the n = 1 cut 
(blue) which means that the configuration is very close to the line Ci. 



Two-Cut Solution, n = 1,3. Two-cut solutions with non-consecutive mode numbers 
allow for even richer properties. As is explained in Sec. 6.1, for some values of filling 
fractions both cuts can develop condensate cores before intersecting with each other. At 
higher fillings the cut with the higher mode number, which initially lies further left, can 
pass through the condensate core of the other cut. As explained in Sec. 6.1, this should 
not affect the inner structure of the cuts. 

From the distributions of numerical roots for different values of the fillings in Fig. 48', 
we see indeed that the group of ra = 3 Bethe roots (left, red) reaches the condensate core 
of the n = 1 roots (right, blue) and passes through the condensate.^'' These numerical 
solutions confirm that the passage of a cut complex through a condensate core does not 
affect the structure of the individual complexes, i.e. their shape and their root densities. 

We did not plot cut contours for all configurations in this case, because similarly as in 
the previous example, the configuration undergoes a phase transition. Namely, this series 
of configurations follows the line P/2Tia = 1.5, which is slightly above the critical value 
~ 1.4 so that we are crossing the line C24 of the Fig. 139) When the two fluctuations 



^''in our example, the (red) n = 3 roots replace some of the (blue) roots in the condensate core while 
the latter take the places of the n — 3 roots. This is merely an artefact of the numerical procedure and 
has no impact on the resulting configuration. 
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denoted by the orange diamonds (n = 2, 4) collide, the corresponding spectral curve 



transforms into a three-cut curve, as is expected from the discussion in Sec. |6.1[ Note 
that we have no means to plot the exact contours for three-cut configurations and thus 
these lines are missing from the figures. 

Conclusions. We have presented numerical solutions for various types of configura- 
tions discussed in this paper, and in all cases we observe a perfect match with the 
analytical predictions of the spectral curve. This confirms the validity of the analysis of 
the thermodynamic limit in the previous sections. 

7.3 One-Cut Solution and Finite-Size Effects 

In this paper we mainly focus on the strict thermodynamic limit. However, the numer- 
ical solutions containing a finite number of Bethe roots allow us to easily analyse the 
leading correction to the thermodynamic limit. In string theory they correspond to the 
semi classical energy shift, see 1 30. 311 for a general discussion of canonical quantisation 
for the spectral curve. 

The energy shift was studied analytically for one-cut configurations of the Heisenberg 
spin chain in |20jl23|l22j|. The result reads 



^(1) = 47r2 



J2lf,'Jl-^'^^^^^^f^-e-2nMl-a)]-nMl-c^) 



fc=i 



(7.4) 



Note that this result equals the sum of zero-mode energies of the fluctuation modes in 
(3.13) upon zeta-function regularisation of the sum. The result may in principle include 
exponentially suppressed contributions, see |32|. It is however conceivable as well that 
these corrections are absent in the Heisenberg/Landau-Lifshitz spectral curve which is 
considerably simpler than the spectral curve for the principle chiral model/string theory. 
The above energy shift was obtained assuming that the roots are distributed along a 
single cut and it is a priori questionable if this result remains true for the case a > Ocond 
when the loop cut appears and the distribution of roots is more sophisticated. We can 
however use the intuition gained studying the strict thermodynamic limit: We saw that 



the energy is described by the same analytical expression (3.3) before and after the 
condensate is formed. Therefore we expect (7.4) to be valid also in the full range of a's. 

Here we perform a numerical test of this natural guess. In Tab. 3 we give a series 
of deviations of numerical values of the energy from their asymptotical values (3.3) for 
different numbers of roots and for two fixed filling fractions below and above the point 
of condensation acond ~ 0.21. To compute 6Eoo we fit the values of the energies from the 
table for finite M by a constant and inverse powers of M. The value of the constant is 
precisely SE^o given in the table. For both filling fractions we see a perfect match with 
the value given by (7.4), which shows that the correction (iTA) is also valid for fillings 
beyond Ocond- 

A possible explanation for this could be based on the algebraic curve. The algebraic 
curve appearing in the thermodynamic limit can be modified to reproduce also the leading 
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a 


5E^o 


5E20 


5E30 


5E,o 


5E,o 


SEqo 


5E^ 


EW 


0.10 


2.372 


2.356 


2.351 


2.349 


2.347 


2.346 


2.3411 


2.3411 


0.33 


-2.273 


-2.472 


-2.539 


-2.573 


-2.594 


-2.608 


-2.6768 


-2.6769 



Table 3: Deviations 6 Em of the numerical energy from its asymptotical value 
for different numbers M of Bethe roots and for two fixed filling fractions a above 
and below Ocond ~ 0.21. 5Eoo is the extrapolation to an infinite number of roots. 
The value E'^^' (|7.4|) is given for comparison. 
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Figure 49: Typical configurations of the two main types of two-cut solutions. 
When the fillings ai^2 are small, the two cuts are disjoint (a), When the fillings 
increase, the cuts attract each other and eventually collide, forming a condensate 
at their common centre (b) . 



quantum corrections |22|24j. The same type of analysis we have performed for the leading 
order hopefully can be applied for this corrected curve as well. 



8 Conclusions and Outlook 

In this article we have considered the general solution (Sec. z) for the spectrum of 
the Heisenberg ferromagnet in the thermodynamic limit by means of spectral curves. 
Equivalently, we have studied the semiclassical Landau-Lifshitz model on S"^, see App. A| 
In particular, we have constructed the general two-cut (elliptic) solution (Sec. ^ and 
studied its properties (Sec. |5j[6j). This allowed us to shed light on the issue of quantum 
stability of classical solutions (Sec. 3,5) brought up in |20,22|. We have also tested our 
claims by constructing some numerical solutions (Sec. 7). Let us summarise our findings 
briefly: 

The two-cut solution is specifled by two mode numbers ni,2 and two flllings ai^2- Its 
quasi-momentum is an elliptic function which encodes all the relevant classical observ- 
ables. In order to quantise this solution, one needs to understand the positions of the two 
branch cuts which correspond to the distribution of Bethe roots. We have studied these 
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(a) a < ac 



(b) ttcond < a < ttcrit 



(c) Qfcrit < a 



Figure 50: The three quahtatively different types of one-cut solutions. For 
small fillings, the single cut is manifestly stable (a). At intermediate fillings, 
the cut obtains a condensate and an additional closed loop cut with a cusp; a 
potential instability is hidden in the condensate (b) , When the filling becomes 
large, the classical solution acquires unstable fluctuation modes and the closed 
loop cut has two cusps (c) Nevertheless, the configuration is physical, albeit 



with non-minimal energy; see the discussion in the main text. 



in detail and we have found a rich "phase space" (Fig. 34p9 ) of different configurations 
depending on the fillings. 

There are two main configurations, see Fig. 49: When the fillings are sufficiently small, 



there will be two separate branch cuts, see Fig. 49a When the fillings are increased, the 
branch cuts bend towards each other so that they eventually would overlap. Instead of 
overlapping, a straight vertical condensate of Bethe roots with constant density forms. 
The configuration looks like a Bethe string but with four curved tails of lower density, 
see Fig. 49b, 

Classically there is no quantitative difference between the two types of solutions, but 
the condensate serves several purposes for the quantum model: It sets an upper bound 
to the density of Bethe roots which is an important feature for quantum stability. In 
some sense, the condensate "hides" some potentially dangerous distributions of Bethe 
roots on an unphysical sheet of the spectral curve. Moreover, solutions with different 
mode numbers (but equal rii —722) can be connected analytically via the condensate. For 
example, there is only a single connected phase space for solutions with consecutive mode 
numbers (n2 = ni -|- 1). Finally, we have observed that the naturally continued overall 
momentum P never crosses the boundary of the first Brillouin zone; it is apparently 
bounded by \P\ < n and the condensate formation plays an important role. At P = ±7r 
the two-cut solution (with consecutive mode numbers) flattens out onto the imaginary 
axis where it connects smoothly the two regimes with P > and P < 0. 

Concerning quantum stability we find that an instability does not actually put a 
threat to the quantisation of some classical solution; the instability is rather caused by 
an unsuitable way of looking at the classical solution. This can be observed conveniently 
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(a) a < ac 



(b) ttcond < a < ttcrit 



(c) Qfcrit < a 



Figure 51: The configurations in Fig. 50b 50c should be understood as degen- 
erate cases of two- and three-cut solutions. This can be seen by the addition of 
small but macroscopic numbers of Bethe roots to the fluctuation points, which 
"regularises" the configurations. The resulting configurations in (b) , (c) are per- 
fectly well-behaved two- and three-cut solutions. Exciting the fluctuation points 



of the solution in Fig. 50c lowers the energy, which is consistent with the (only 
seemingly) unstable excitation modes. 



by means of the one-cut solution (Fig. '50^ Here there are three main classes of cut 
configurations (which serve as subclasses for the two-cut solution): 

For small filling, there is just a simple cut, see Fig. 50a, in which case the solution is 
perfectly stable. For intermediate filling the naive one-cut solution violates the stability 
criterion but without any observable flaws like unstable fluctuation modes. In this case, 
the physical cut conflguration develops a condensate and a loop with a cusp, see Fig. 50b| 
which achieves to bypass the stability criterion; it is nevertheless semiclassically equiv- 
alent to the one-cut solution. For large flUing, the naive one-cut solution acquires an 
unstable fluctuation mode which appears to render the solution unphysical. The physical 



cut conflguration now has a condensate and a loop with two cusps, see Fig. 50c which 
gives a clue of how to understand the instability: 

The latter two dressed conflgurations should be interpreted as singular limits of two- 
cut or three-cut solutions, respectively, see Fig. 51, In the naive one-cut conflguration, a 
fluctuation would break apart the loop at a cusp by a tiny amount. This is flne when the 
cusp is on the real axis (for intermediate fllling), but away from the real axis (for large 
flUing) it leads to complex momenta and thus an instability. Conversely, fluctuations of 
the regularised three-cut conflguration would merely change the fllling of either the two- 
cut complex or the third cut; in either case the conflguration remains symmetrical with 
respect to the real axis, and there is no physical instability. Therefore the classical one- 
cut state with large fllling does have a corresponding quantum state. This state however 
does not minimise the energy locally in accordance with the apparent instability. The 
energy is lowered by taking away Bethe roots from the third cut and moving them onto 
the two-cut complex. At the end of this process the third cut has evaporated and only 



the two-cut complex remains which is a genuine two-cut solution (Fig. |52c|l0e). 
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(a) a < ac 



(b) ttcond < a < ttcrit 



(c) Qfcrit < a 



Figure 52: The proper physical continuation of the one-cut solution beyond 
the critical filling acrit is the two-cut solution, Fig. |10e[ Unlike the degenerate 
three-cut solution, Fig. 50c, this configuration locally minimises the energy. The 
phase transition from the one-cut to the two-cut solution occurs when the right 
fluctuation point collides with the cusp of the closed loop cut. Shown here is the 
regularised version of this process, where the fluctuation points carry a small but 
finite number of excitations. 



This leads to the conclusion that the one-cut solution undergoes a phase transition 
at the point which separates intermediate and large filling: When the one-cut solution 
is analytically continued from small to large filling one obtains a degenerate three-cut 
solution. However in the vicinity of the latter there exists a two-cut solution with slightly 
smaller energy. It is useful to consider this solution the proper physical continuation of 
the one-cut curve. We find that at the phase transition point the second derivative of 
the energy is discontinuous. An exact expression for the discontinuity has been found 
on the basis of numerical data; an analytic confirmation of its form would be desirable. 
Technically, the phase transition turns out to be equivalent to the one observed in 2D 
QCD on S"^ 1 29 1 or equivalently a particular matrix model. In App. B we summarise 
the relevant functions and illustrate the phase transition for that model. Note that the 
apparent instability and the phase transition are merely artefacts of the thermodynamic 
limit in accordance with general results of statistical physics; they are resolved naturally 
at finite size, i.e. by the quantum theory: At finite L there exists no precise notion of 
continuous cuts and the distribution of Bethe roots onto these cuts is often ambiguous. 
Thus two cuts never actually end on the very same point which was an essential assump- 
tion for the instability and the phase transition. In particular we show that numerically 
exact solutions with finite L are smoothly interpolating between two phases and there is 
no real instability when the condensate forms. 

We hope to have presented a fairly complete picture of the Heisenberg ferromagnet 
in the thermodynamic limit, its quantum fluctuations and the issue of stability. However 
there are several further questions which may be addressed in future work. One of them 
would be how precisely the third cut interacts with the two-cut complex in solutions of 
the type shown in Fig. 50cg51c, An investigation of this issue would require an explicit 
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expression for the (hyperelliptic) three-cut quasi-momentum. Further, it would be in- 
teresting to construct the general elliptic solution to the Landau-Lifshitz model which 
includes the solutions obtained in |15|33|. It should be equivalent to the two-cut solution 
presented here and the circular solution (A. 8) should transform into the elliptic solution 
in the same way the two-cut solutions emerges from the one-cut solution at large fillings. 
Another question would be how the thermodynamic solutions of the Heisenberg magnet 
behave under the influence of a magnetic field and how this would change the picture of 
the phase space. 

Concerning the AdS/CFT correspondence, a generalisation of the results to the SU(2) 
principal chiral model describing the dual string theory on M x S^ would be desirable. 
As discussed in |18|, the solution to this model is given by a spectral curve that reduces 
to the Heisenberg spectral curve in a certain limit. Because of their strong structural 
similarity, we expect that stability and loop-formation work in qualitatively the same 
fashion. The phase space of general two-cut solutions (or elliptic string states) is likely 
to be similar as well, but new features may appear due to the more elaborate singularity 
structure on the curve. Therefore it is worth obtaining and discussing the general two- 
cut curve for the principal chiral model. Some special cases of two-cut solutions have 



appeared in 27,34 . Analogously one could find the corresponding elliptic solution of 
the equations of motion explicitly, e.g. using the results of |30|. This solution may be 
related to the solutions studied in |35|. 

Finally, it remains an open question whether a similar study of the thermodynamic 
phase space can be carried out for spin chains which have different symmetry groups, 
quantum-deformed symmetry, or for which the spins transform in different representa- 
tions. For example, studying the spectrum for symmetry groups of higher rank requires 
working with spectral curves with more than two sheets, which complicates the analysis 
substantially. Moreover for more general symmetry groups the algebraic curves are quite 
sensitive to different boundary conditions of the spin chain |24|, thus making the analysis 
even more complicated (and interesting). 
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A The Landau-Lifshitz Model 

The long-wavelength limit of the Heisenberg ferromagnet is equivalent to the Landau- 
Lifshitz model on S^ |16,15|, see also |36| for a review of Landau-Lifshitz models in 
connection with the AdS/CFT correspondence. Here we shall discuss the simplest cir- 



68 



cular solution and its fluctuation spectrum which corresponds to the one-cut solutions 
discussed in Sec. El 

A.l Definitions 

The Landau-Lifshitz model on S"^ is a classical two-dimensional sigma model with two 
fields 6{a, r) and 0((t, r) which are the standard angular coordinates on S"^. We identify 
the site k of the corresponding Heisenberg chain with the coordinate a = 2'n'k/L of the 
sigma model. The expectation value of the spin projection Sz at that site corresponds 
to the value | cos 9 at a. 

In analogy to the closed spin chain we take the world sheet to be closed, i.e. the 
coordinate a is identified with a + 27r: Consequently, the embedding coordinates 6 and 
(f) (more precisely the equivalence classes describing points on S"^) must be 27r-periodic 
in a. This leads to the following periodicity conditions 

e{a + 27T,T) e±e{cx,T) + 27TZ, 0((T + 27r,r) e 0(a,r) + 27rZ+^7r(l^l). (A.l) 

The Lagrangian reads 

C = -^^cose-^(e'^ + (j)'^smH), (A.2) 

47r 2L 

where the constants L/4n and 7r/2L are adjusted to match with our conventions for the 
Heisenberg chain. In fact the first term in the Lagrangian originates from a topological 
Wess-Zumino term in three dimensions and for a fully consistent treatment one should 
make the replacement 

1 a f^f d(f)dcose d0dcos^\ 

The equations of motion following from the Lagrangian are 

0= f'^j (26' (f)' cos 6 + 0" sin 6) , 

<P=('^) {(P'^ cose -9" CSC 9). (A.4) 

Taking the derivative of C with respect to 9 and (p yields the conjugate momentum 
densities 

Tig = , 11(1, = cos6'. (A. 5) 

An 

The Lagrangian is homogeneous and therefore the total energy and total momentum 

are conserved. The rescaled total energy and a suitably shifted total momentum read 

E = EL = L I {■Ke9 + -K4>^-C)d(j= | f [9'^ + ^'^ sin^ 9)da, 

P=^f^'da+'^f {ng9' + 7r^0') da = ^ f {l - cos9)(P' da. (A.6) 
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Furthermore the Lagrangian is invariant under shifts i— * (j)+e which yields the conserved 
total spin S* = L(| — a) with the filling 



2 + ;^ y 7r<^da = — y (1 - cos^) da. 



a = - + - n^da = —- (l-cos9)da. (A. 7) 



A. 2 Circular Solution 

The simplest non-trivial solution of the Landau-Lifshitz model is the circular string at 
constant latitude 6q with n windings |'18'| 

e{a, t) = 00, (f){a, t) =n(T + (27r/L) Vr cos ^o- (A.8) 

It has the same charges as the one-cut solution of the Heisenberg magnet ( |3.3[ ) 

a= Kl-cos^o), P = 2Tma, E = Air'^n^ail - a). (A.9) 

The momentum P and the filling a are defined such that they both vanish for the 
vacuum solution at 6^0 = 

a = P = E = 0. (A.IO) 

A. 3 Fluctuations 

We would like to obtain the spectrum of fluctuations around this solution. For the case 
of a = ^ this problem was solved in [37] which coincides with the ultra-relativistic limit 
of spinning strings discussed in |21|. To this end, fluctuation modes with Fourier mode 
k and frequency Uk are added to the background solution with the following ansatz 

69 = —^ cos(^ka + UkT + c) , 6(f) = —^ sm(^ka + UkT + d). (A. 11) 



From the equations of motion it follows immediately that (i = c as well as 

/2vr\V , . ^,2 .\ ^27r\V , , aP-n^sin^^ 
cu^ = -- Ink cos 6' + - fc^ sm 6* = -- Ink cos Q + 



L J \ a J \ L J \ b sin 6* 

(A.12) 
This flxes the frequency and ratio of amplitudes 



,'27iY / , , ^ ,9 / 4n2a(l -a) 

^k= [^] { 2nk{l - 2a) + k\ I ^ ' 



L \ ^ 'V A;2 



^"'"('-"l (A.13) 



.2 



a 2^a{l-a) V F 

Now it turns out that the total energy, total momentum and flUing are not yet 
completely flxed. They still depend on the value of a and curiously also on the second 
variation of the a- independent mode^ 

6^6 = {. (A. 14) 

ij 



^^This feature is related to the fact that this mode is already macroscopically filled in the classical 
solution. Essentially / corresponds to changing a in the classical solution by a microscopic amount. 
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We shall parametrise our ignorance by two constants Nk and Nq as follows 



a 



V2N~J {l 



aI Aii?a{l — a) 



^ 4/^ _ 4n2a(l - a) 



2Ja{l-a) V P 



This leads to the following shift in the observables 

.2 ^0 + Nk 

^^P = iVo^ + iV,?!^, (A.16) 



,,~_,Ar:^n\l-2a) , ^ 4vr2 / ^^ , ^,^^^^ ^^^ , . /^ 4^^! - «) 



6'E = No \ '- + Nk—- \n{n + 2k){l-2a) + k\l 



The excitations enumerated by Nq and A'^^ carry each one quantum 1/L of filling equiv- 
alent to one unit of spin. They have mode numbers n and n + k, respectively.^^ Their 
energies can be read off from the above expression. It is not hard to see that a quantum 
of mode n essentially just increases a by 1/L in the classical solution, cf. (A.9). The 
charges for mode n + k are in complete agreement with the fluctuations in the ferromag- 
netic Heisenberg magnet (3.13). A careful semiclassical analysis actually shows that the 
numbers A'o and Nk must be integers. ^^ 

The classical solution (A. 8) is unstable if some of the oscillator frequencies Uk are 
complex. Obviously, this is the case if 

2n^/a{l-a) > 1. (A.17) 

The transition from stable to unstable classical solutions exactly happens at the critical 
density Ocrit (3.15), where the one-cut solution of the Heisenberg magnet transforms into 
a two-cut solution, as discussed in Sec. ,5.4, 

B The Douglas— Kazakov Transition 

Here we shall review the Douglas-Kazakov transition |29| for QCD on S'^ on a technical 
level. It mimics the transition between our one-cut and two-cut solutions, albeit in a 
somewhat simplified and more transparent fashion. 



^^The shift by n w.r.t. the Fourier mode in 69,S(j) is induced by the linear dependence of (f) on a. 

^•^Note that our Lagrange function evaluated on the classical solution is not zero and does depend on 
the filling a. Unfortunately, this obscures the correct phase space normalisation in connection with the 
shift 6^9. If one uses the Wess-Zumino term (A. 3) instead, the classical Lagrange function is zero and 
the above problem should be absent. 
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(a) rational regime, A < Adit (b) elliptic regime, A > ^crit 

Figure 53: Density function p(/i) in the respective regimes. 

B.l Definitions 

The description of the partition function in the large- A^ hmit of QCD on S"^ uses a density 
function p{h) with support on the real axis. It must obey certain integral equations which 
will not be needed in detail here. 

In order to extract the relevant observables it is convenient to introduce the contin- 
uation u{h) of p{h) into the complex plane. The function u{h) has a branch cut on the 
real axis such that the real/principal part equals the density 

p{h) = lu{h + ie) + lu{h-ie). (B.l) 

The function u{h) is the analog of our quasi-momentum p{x) and its derivative u'{h) is 
an algebraic curve. Apart from the branch cut, the only singularity is at /i = cxd with 



u{h) 



^(-4" + | + ^(f' + i) + 0(l/"=)) 



(B.2) 



This expansion also defines the area A and the derivative of the free energy F'. It is clear 
that these quantities play similar roles as the total filling a, momentum P and energy 

E. 



B.2 Rational Regime 

The simplest solution has one branch cut at [— /^max, +^r 
quently is of semi-circle form 

u{h) = ^^A{4-Ah^), 
with /imax = 2/yZ, see Fig. |53a[ 



The function u{h) conse- 
(B.3) 



The resulting free energy from (jB.2|) is simply [38 

1 1 

2A~ 24' 



F' 



(B.4) 



This function is physically correct for A k. {}. However, it was noticed that it does 
not have the correct behaviour at A ^ oo. For instance -^' ^ ^^ but not F' -^ 0. 
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B.3 Elliptic Regime 

On the technical level, the pathology can be attributed to a density which grows beyond 
a threshold of 

Pmax = 1- (B.5) 

A density with p{h) > 1 is unphysical. Indeed for A > Acrit with 

^crit = vr^ (B.6) 

the maximum density p(0) = VA/n exceeds the maximum allowed density p(0) > pmax- 
This problem is resolved by inserting a condensate with constant density p{h) = 1 on 
the interval [—&,&] |29|. The density then falls off to zero on the intervals [—a, —b] and 
[b, a] . This leads to the two-cut solution, cf . Fig. 53b| 

u(h) = -^J(a^-h^W-b^) U(byh^\ b^a^). (B.7) 

irah 

The expansion ([372) at h ^ oo leads to 



(B.8) 



where 



A = 4K(g)(2E(g)-(l-g)K(g)), 
p; ^ 4(l + g)E(g)-(l + 3g)(l-g)K(g) _ ^ 
24(2E(g)-(l-g)K(g))' 24 

1 ._ v^ 



(B.9) 



2 E(g) - (1 - g) K(g) ' 2 E(g) - (1 - g) K(g) ' 

As usual this implicitly defines F' as a function of A. 
Now for A ^ oo one finds 

F' = e-^/2 + ( 1 ^2 - 3A+l)e~^ + (|^^ - f ^' + 14^^ - 8A + 4)6-=^^/^ ^ 0(e-2A) (g^^Q^ 
in agreement with physical expectations |39|. 

B.4 Transition 

We can now compare the two regimes. The expansion of the derivative of the free energy 
around the critical point A = Acrit = tt^ yields in the two cases 

p/ _ 1 1 A — Acrit _|_ (^ ^ ^crit)" [A — Acrit) 

■"^^ ~ 2^ ~ 24 27r4 ^ 2^e ^8 + • • • ' 

p, _ 1 1 A-Arit MA-Aritf IHA-Aritf .^^^. 

^^^^"27r2 24 27r4 + 2vr6 4n^ +•••' ^^'^^^ 

see Fig. 54a. The deviation is therefore of second order 
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(a) F'iA) 



(b) Comparison of p{h) 



Figure 54: Transition of the derivative of the free energy F' at A = ^crit = t^"^ 
and comparison of density functions shghtly above ^crit- 

implying a third-order transition in the free energy F. 

In fact the function u{h) deviates only quadratically at a generic point away from the 
branch points, see Fig. 54b, It can be expected that the same holds whenever an elliptic 
curve degenerated into a rational curve. This leads to the conclusion that all quantities 
which can be read off directly from the curve (e.g. in an expansion) will deviate at least 
quadratically. 



C Integration of the Partial Filling Fractions 



Analytically integrating the partial filling fractions ai, 02 of the two-cut solution (4.11) is 
possible because they can be expressed as integrals over A-cycles ( ,2.23, ) of the derivative 
Po(x), which is an algebraic function: 



1 



"1,2 



2Tri 

1 
27n 



2ni 



1 



2m 
1 

2TTi 



p{x) dx 

Po o fJ'{x) + C dx 
iiPo + C) o fi){x)dx 
{po + C){y){f,-'y{y)dy 

M(^1.2) 

PQ{x)fL~^{x) dx 



Al,: 



Ai 



I alibi -x^ 



27r« Mai,2) V ^o(«o - ■^^) 
An {bl{sx — uY E(g) — (ogS — ux){hls — ux) K(g)) 
aoirx — t){sx — u){hl — x"^) 



dx . 



(C.l) 
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In the following, this integral will be evaluated for the case < ao < bo < oo. The 
result can then be analytically continued to the case of complex oq and bo = oq. For 
< ao < &o < oo the contours Ai^2 encircle the real intervals Ji := [ao,&o], h '■ = 
[—bo,—ao]- Except for an overall sign, the integrand Pq{x)^~^{x) has the same values 
above and below the intervals /i_2- Therefore, the integrals along the contours Ai^2 can 
be expressed as integrals over Ii^2'- 

f Po{x)fi^^{x)dx = 2\im. / pQ{x)fi^^{x)dx . (C.2) 

m(^i,2) h,2-i£ 

Since 



fog - (x - ie)^ \ £^0 J < for X G /i 
ag - (x - iey J I > for x G /2 



im( i r ^.^(j ^r: r^^v, (c.3) 



Iq — (x — iey 
the square root in the integrand on the right hand side becomes 



bl-{x-iey e^o^f-i (x G /i)\ j bj - x^ 



ag - (x - ie)2 \ i (x E h) j y x"^ - a^ 



(C.4) 



The partial filling fraction integrals can now be written as 



' ^boJj,^,^{x^-al){bl-x^)\ x-l X-"' ' ^ ' 

where ai has the positive and 0:2 the negative sign. The constants are 

^^_ u'K{q)-bls^E{q) 

rs 
{airs — tu){blrs — tu) K(g) — bl{st — r-u)^ E(g) 



At/r 



r'^{st — ru) 
{ajs^ - u''){bls^ - u')K{q) 
s'^{st — ru) 



The integral (C.5) can be performed term by term. The first term is the simplest: 



^ dx = -^ K(g) . (C.7) 

j^^^ ./(^^^^^^)(^[^^) bo 

The second and third term are of the same form, but the results for either contour differ 
by a sign: 



dx 



Ji,^ix^-al)ibl-x^)x-c 

X + c 



dx 

ao V(x2-ag)(6g-x2)(x2 -c2) 
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TT 



c^-a" 



al ~ 



2\/(«o-c2)(6g-c2) bo{c^-al)\ c^ 

dx 



K(g) + i n(c 



/2 



v/(x2-a2)(62_3;2)x 



TT 



c^-a^ 



an ~ 



V(«o-c2)(&o-c2) &o(c2-a2)V c2 
For obtaining these two forms, the integrals 

rbo 



K(g) + ^n(c 



(C.8) 



X 



dx 



ao 



y/{x^-alM-x^){x^-c^) Jal 2^{t~al){hl-t){t-c^) 



dt 



Tl 



2^{al-c^){hl-c^) 



(C.9) 



and 



60 



ao 



\/{x'^ — ao)(feo ~ x'^){x'^ — c^) 



dx 



60 (c2 - ag) 



c^-a" 



K(g) 



^nfc^ 



_ (CIO) 

were used. The latter can be found in |40|, formula 217.12. Putting it all together and 
writing the occurring square roots in the form (4.17) in order to select the branch of the 
square root consistently, one arrives at the result ( j4.19, ). 

D Integral of the Two-Cut Quasi-Momentum 

In order to find the physical cut contours for a given two-cut quasi-momentum p(x), 
one can use the fact that expression (4.22) must be real on the physical contours. That 
expression contains the integral of the quasi-momentum, which will be given here. Direct 
integration yields 

iaoAn 



-1 hl{hl-y^ 






rs{hl - y2) y a2 y hl{al - y^) 

■ {bls^ E(g) - u^ K(g)) F | i arsinh 



/)2 




1-q 



ao{bo + y) 



a-oibo - y) 



+ 



4Araao V ("0 - &o)(y - ^o) V (^o + &o)(ao - v) 



ft2 



(ao-6o)(ao + y) al{bl - y"^) 



(ao + bo){ao -y)]] b^a^ - y'^) 

bl{st — ru^ E(g) — {airs — tu){blrs — tu) K(g) 
(ttgr^ — t'^){st — ru) 



^vAy) 
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st — ru 

bls{st — ru) E(g) + (aorsti — tu"^ + aQr{u^ — hls^)) K(g) 
2aors(aor — t) 



• i:* I arcsm \ / -— ^ 



(oo - &o)^ 



(D.l^ 




(ao + bo) {aoA - B) . / (ao - feo) (ao + a^) 

-— ; arcsm -i / -— 

(ao - Oo)(ao^ + B) V (ao + Oo)(ao - x) 



(ao + bo) 



{ao - 60)^ 



(D.2) 



and where F[z\q) is the incomplete eUiptic integral of the first kind: 



/■sin 2 -1 

F{z\q)= ^dt. (D.3) 

Jo ^{l-t^){l-qt^) 

Using the integral (D.l), one can obtain a closed expression for the integral A(x) of the 
two-cut quasi-momentum p{x) ( |4.11| ), which can be expressed as 

/X PX 

Viy) dy = xp{x) - j yp\y) dy 

/x ^iiix) 

yPoiKy))f^'{y) ^y = xp{x) - / fi'^{z)pQ{z) dz . (D.4) 
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